Control of dipolar relaxation in external fields 
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We study dipolar relaxation in both ultra-cold thermal and Bose-condensed chromium atom gases. 
We show three different ways to control dipolar relaxation, making use of either a static magnetic 
field, an oscillatory magnetic field, or an optical lattice to reduce the dimensionality of the gas 
from 3D to 2D. Although dipolar relaxation generally increases as a function of a static magnetic 
field intensity, we find a range of non-zero magnetic field intensities where dipolar relaxation is 
strongly reduced. We use this resonant reduction to accurately determine the S = 6 scattering 
length of chromium atoms: a§ = 103 ± 4ao- We compare this new measurement to another new 
determination of a§ , which we perform by analysing the precise spectroscopy of a Feshbach resonance 
in <i-wave collisions, yielding ae — 102.5 ± QAao. These two measurements provide by far the 
most precise determination of a§ to date. We then show that, although dipolar interactions are 
long-range interactions, dipolar relaxation only involves the incoming partial wave I = for large 
enough magnetic field intensities, which has interesting consequences on the stability of dipolar Fermi 
gases. We then study ultra-cold chromium gases in a ID optical lattice resulting in a collection of 
independent 2D gases. We show that dipolar relaxation is modified when the atoms collide in reduced 
dimensionality at low magnetic field intensities, and that the corresponding dipolar relaxation rate 
parameter is reduced by a factor up to 7 compared to the 3D case. Finally, we study dipolar 
relaxation in presence of radio-frequency (rf) oscillating magnetic fields, and we show that both the 
output channel energy and the transition amplitude can be controlled by means of rf frequency and 
Rabi frequency. 

PACS numbers: 34.50.-s, 67.85.-d, 37.10.Jk 



Up to recently, only short-range isotropic interactions 
have played a significant role in the physics of ground 
state cold atoms. In contrast, dipolc-dipole interactions 
are long-range, and anisotropic. Strong dipole-dipole in- 
teractions have profound consequences on the properties 
of either ultra-cold thermal or quantum degenerate Bose 
or Fermi gases. Dipole-dipole interactions greatly modify 
the stability properties [l( and the excitation properties 
Q of Bose condensed gases, and can even lead to new 
quantum phases H, or novel types of superfluidity 
One may also use the long range character of dipolar 
interactions to entangle two [5| or more atoms, with in- 
teresting perspectives for quantum computation Q. 

Strong dipole-dipole interactions arise when an atomic 
or molecular species carries a strong permanent magnetic 
or electric dipole moment. Good candidates therefore in- 
clude heteronuclear molecules with large electric dipole 
moments (which were recently produced at large phase- 
space densities 13), or atoms with large electronic spin 
(so far erbium Q, dysprosium Q and chromium). Up 
to now, chromium is the only species with large dipole 
moment for which a quantum degenerate gas has been 
produced 1(J U|. Smaller dipolar effects were also ob- 
served in a BEC of potassium for which the scattering 
length can be precisely tuned to zero by means of a Fes- 
hbach resonance (l2| . 

Interesting new physics comes at play when one also 
considers the spin degree of freedom. Spin dynamics of 



optically trapped multi-component Bose-Einstein Con- 
densates (also known as spinor condensates) has been 
observed Coherent oscillations between the spin 

components is driven by centrally symmetric short range 
exchange interactions, and the total magnetization in the 
system is conserved. In [l5[ a first dipolar effect was ob- 
served on the spin texture of a Rb spinor condensate. 
Dipole-dipole interaction will introduce additional new 
features in spin dynamics as it couples the spin degree of 
freedom to orbital momentum, due to its asymmetry. For 
example, when the magnetic dipole mean-field created by 
the atoms locally dominates the applied static field, spin 
precession occurs and the initial magnetization of the sys- 
tem is transferred into mechanical rotation (the equiva- 
lent of the Einstein-de-Haas effect), with appearance of 
structures similar to vortices [16j |. Spontaneous circula- 
tion was also predicted for some experimental parameters 
of ground-state spinor dipolar Bose-Einstein Condensates 
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Coupling between spin and angular momenta is pro- 
vided by dipolar relaxation, a two-body process in which 
the total magnetization changes. The main topic of this 
paper is to discuss how to control both the dipolar relax- 
ation rate in cold dipolar gases and the energy released in 
a single dipolar relaxation event. Controlling (and pos- 
sibly reducing) the rate of dipolar relaxation is interest- 
ing, as this inelastic process generally limits the lifetime 
of dipolar gases when all particles are not in the low- 
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est energy state, and in particular for spinor condensates 
studies. (Similar limitations occur for ultra-cold Rydberg 



atoms |18j). In addition, the energy released in a dipo- 
lar relaxation event is set by the Larmor frequency of the 
atoms, typically much larger than the chemical potential. 
To observe coherent spin dynamics due to dipolar relax- 
ation in spinor condensates, the Larmor frequency should 
be of the order of the chemical potential or less. As it 
is technically difficult to control weak magnetic fields at 
the required precision, it is appealing to device new tools 
to control the energy in the output channel of dipolar 
relaxation. 

In the first part of this paper, we uncover a range of 
magnetic field strengths for which the dipolar relaxation 
rate parameter of S — 3 Cr atoms in the stretched high 
field seeking state mj = 3 is strongly reduced, both for 
a thermal gas and for a BEC. This reduction stems from 
the last node of the s— wave collision S — 6 radial wave- 
function. We analyse the resonant reduction of dipolar 
relaxation as a function of magnetic field, which leads 
to a new determination of the S = 6 scattering length 
of Cr ag. Due to the proximity of a shape resonance in 
S = 4, 1 = 2, the value of the S = 4 scattering length 0,4 
also has a significant contribution to the dipolar relax- 
ation rate, and our analysis therefore also leads to a new 
estimate of 04. 

Then, we demonstrate that dipolar relaxation in both 
a BEC and a thermal gas is only due to s-wave collisions, 
provided the magnetic field is large enough. This is in 
contrast to elastic dipole-dipole interactions, to which all 
partial waves contribute due to the long range character 
of dipolar interactions. 

We also describe experiments at low magnetic fields, 
where we use optical lattices to strongly confine the mo- 
tion in one direction of space. In the 2D cold gases that 
we thus produce, the density of states at the energy of 
the output channel is strongly reduced compared to the 
3D case, which thereby leads to a strong reduction of 
dipolar relaxation in the ms = 3 state. In addition, 
dipolar relaxation in the ID lattice depends on the ori- 
entation of the magnetic field relative to the plane of 
the 2D direct consequence of the interplay be- 

tween the anisotropy of the trapping potential and the 
anisotropy of dipole-dipole interactions. 

Finally we describe experiments where dipolar relax- 
ation of atoms in a BEC of 7715 = — 3 Cr atoms (usually 
energetically forbidden) is triggered by an rf photon. We 
show that the rf-assisted loss mechanism is due to dipolar 
relaxation between rf manifolds, and that the loss param- 
eter and the energy released in a dipolar relaxation event 
can be controlled by means of the rf Rabi frequency Q 
and the rf frequency oj. 



THE EFFECT OF MOLECULAR POTENTIALS 
ON DIPOLAR RELAXATION 

Dipolar relaxation has already been experimentally 
studied in Li [ljj, metastable He [jof . and Cr gases. 
Dipolar relaxation cross sections have also been theoret- 
ically estimated by means of coupled channels methods 
22], or using first order perturbation theory In 
chromium, dipolar relaxation has been estimated using 
the bare first order Born approximation ignoring all ef- 
fects of the short-range molecul ar p otentials, and com- 
pared to experimental results in [21]. Here, we will show 
that while this approach is valid at relatively low mag- 
netic fields, it is crucial to take into account the molecular 
potentials at higher magnetic fields. 

In this part, we study the specific example of two 
\m,g — 3) colliding chromium atoms, in a magnetic field 
of intensity B, set along the z direction. S — 3 is the 
spin of ground state chromium atoms. Dipolar interac- 
tions between two particles separated by r can be written 
as: 
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where d 2 — Ho(gsfJ-B) /47r and r = \\r\\. fiQ is 
the magnetic constant, g$ ~ 2 is the Lande factor 
for ground state chromium atoms and [Ib is the Bohr 
magneton. Starting with a pair of atoms in state 
\ms = 3,m5 = 3} = |0), there are two relaxation chan- 
nels due to the dipole interaction operator Vdd, channel 
1: 
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and channel 2: 
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A.EU=(i> 2 )) are ^ ne g ams m kinetic energy by the pair 
of atoms after a dipolar relaxation event for both chan- 
nels. For atoms in the state mg — 3, dipolar relaxation 
is exo-energetic: the kinetic energy increase is directly 
proportional to the Larmor frequency. In contrast, dipo- 
lar relaxation of chromium atoms in the lowest state of 
energy 777,5 = — 3 is endoenergetic, which has been used 
to cool dipolar atoms 24 1. 

If we characterize the relative motion of the particles 
by the relative wavefunctions ^i n (r) and &out(r) for the 
incoming state and the ouput state respectively, the two 
channels are described by the following matrix elements: 
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x, y, z are the relative coordinates of the two atoms sep- 
arated by f. V 7= (i i 2) is the matrix element for dipolar 
relaxation in channel 1 or 2. As Vi and V% involve terms 
in (a; + iy) dipolar relaxation induces transitions between 
the different partial waves. 

To estimate the cross-sections for dipolar collisions in 
the framework of the first order Born approximation, one 
calculates the Fourier transform of the dipole-dipole cou- 
pling independently for each channel, V(j) (kf — kij , for 

any value of the initial and final wave-vectors ki and 
kf. Then, the total cross-section is obtained by sum- 
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over all possible values of kf and 



averaging over all possible ki. 

The results of the calculation for the elastic and inelas- 
tic cross-sections (for channel 1 and 2) are respectively 
2l|: 
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where 



f(u) = l + e 
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e = 1 for bosons (as in [2l|), and e = — 1 for fermions. ki 

(1) 

and kf are related through the conservation of energy 
equation: 



h 2 (k 
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where m is the atom mass. 

Eq. (j8]) describes the elastic cross section due to dipole 
dipole interactions. It is worth emphasizing that this 
cross section does not vanish at small collision energy for 
fermions, showing that non zero partial waves contribute 
to the elastic cross section even at zero collision energy 
(as polarized fermions do not interact in s— wave), a di- 
rect consequence of the long range character of dipole- 
dipole interactions. In addition, it is enlightening to em- 
phasize that ^/oq is, to within numerical factors, equal to 
the range of the dipole-dipole interaction potential Rdd, 
denned by ^ = 

Eq. © and (TlU|) give respectively the cross-sections 
for dipolar relaxation in channels 1 and 2. As we shall 
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FIG. 1: Dipolar relaxation from the molecular physics point of view, 
illustrated for channel 2 (sketch); the molecular curves of the initial 
/ — potential and the final 1 — 2 potential cross at interparticlc 
distance _R* . 



later see, the experiment does not directly provide a mea- 
surement of cross sections, but rather a measurement of 
the corresponding loss rate parameters. For one particle 
moving at velocity v% in a still medium of density n, the 
loss rate is Tj — 2najVi — 2najh^j^. The factor 2 is 
due to the fact that two atoms are involved: both are 
lost in a single dipolar relaxation event when the gain in 
energy of each atom is higher than the trap depth. To 
relate this rate parameter to the rate parameter in a gas 
of N particles, one should multiply this result by N/2 to 
properly account for all pairs of particle. The final loss 
rate parameter is therefore (3^ = ajm = 2a jU^- 

The main weakness of the calculation shown above is 
that it neglects all other interatomic potentials (i.e. the 
electrostatic couplings). At low collision energies, cen- 



trifugal barriers h 1(1 + l)/mr 2 for particles colliding in 
partial waves I > prevent the particles from getting 
close to one another, so that one can indeed safely neglect 
the electrostatic potentials (Cq/t 6 , Cg/r 8 , ...). However, 
particles colliding in s— wave can approach each other 
and experience the short range interaction potential. To 
take this into account, it is useful to develop the collision 
wave function in the partial wave basis. At long dis- 
tances (i.e. further than the van der Waals range of the 
potential R vd w, defined by h 2 /(mRl dW ) = C 6 /R® dW ), 
there are analytical solutions to the Schrodinger equation 
describing the relative motion of two particles, for any 
partial wave. In a first approach, we will use these ana- 
lytical wavefunctions to calculate the dipolar relaxation 
rate parameter. As we shall see, this approach (although 
inaccurate) will provide some physical insight into the 
mechanism at play in dipolar relaxation. 

For a spherical gaussian BEC, the situation is rela- 
tively simple, as the symmetry of the order parameter in- 
sures that the reduced two-body wave function has a pure 
s— wave character. In a BEC, dipolar relaxation couples 
the purely s— wave collision wave function to I = 2. In- 
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deed, the matrix element for both channels of dipolar 
relaxation are proportional to a spherical harmonics Y£ 
(see eq. (O and Q), and thus couple I — to I = 2 by 
the following angular matrix elements: 

V^\r) = (0 > Z=0|V <w |l,Z = 2,m, = l) 
V 15 r A 

V^ 2 (r) = (Q,l = Q\V dd %l = 2,rr H =2) 



= 2>Sd 2 



2 1 



15 r 3 



(13) 



Having thus determined the matrix elements for the 
orbital coupling from I = to I = 2, we now turn to the 
evaluation of the radial coupling of the corresponding 
wave functions. At zero temperature, and in free space, 
the incoming radial wave function at sufficiently large 
distances (i.e. larger than R v dw) can be written as: 



F in {r) 



(14) 



where vo is a normalizing constant (see below). At large 
distances, the ener gy- n ormalized I — 2 output radial 
wavefunction reads [25j: 
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j2 is the second spherical Bessel function, and fc) are the 
final wave vectors for each channel of dipolar relaxation 
h 2 (k ( f j) ) 2 /m = AEtf). The coupling between F in and 
F ut is given by the Fermi golden rule 261 ]: 



7ir (i) = 2tt 



F m (r)F^\(r)V^ 2 (r)r 2 dr 
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(which can be calculated analytically) is the event 
rate associated to the dipolar relaxation in channel (j) 
[26l ]. It is now important to clarify the value of the 
normalization constant vq. This can be done for the 
example of an isotropic harmonic trap with trapping 
frequency lo/2tt. At long distances, the radial wave 
function is essentially set by exp(— r 2 /4a| fo ), where 
ciho = yft/tjmcJ)) whereas at short distances (but not 
too short, as explained above), it is set by F, n (r). As 
a ho » a 6> one can safely approximate the initial 
I = radial wave function in the trap by F i ^ ap (r) = 
Fi n (r) exp(— r 2 /Aa 2 HO ). Using this form for the incom- 
ing wave function, one can now calculate the coupling 
between F^ ap and F out by an equation analogous to eq. 
([To]) . The integral is still analytical, and almost identical 
to the result of eq. (|16p , for our experimental parameters. 
We therefore use Fi n of eq. [TJ] for the calculation of the 
dipolar relaxation coupling, and F^ ap for the determina- 

2 dr = 1, 



tion of vq. Normalizing F i ™ p using 



we find vq = [clho) (27r) 
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To now relate calculated for a pair of atoms, to 
the value of the rate parameter for dipolar relaxation f3 r , 
we recall the definition of f3 r , through the local equa- 
tion = —(3 r n 2 . Integrating this equation over the 
profile corresponding to the density probability of two 
particles in the vibrational ground state of the trap, we 
find that the associated loss rate is 2T — ^72 Pr (r is 
the collision event rate so that 2r is the loss rate), where 

no = 2 (^) 3 ^ 2 , so that T = ^(3 r . As a consequence, we 
obtain: 



0%* = r(% = ^™sW (1 - kf\ir) 
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In the limit where fe, — > (BEC case) and ae = (no 
molecular potential), these equations exactly correspond 
to half the value of the results of the rate parameter as- 
sociated to the cross-sections given in eqs. ((9|) and (fT0|) 
(corresponding to 2l|). The difference of a factor of two 
arises from the different way of symmetrizing a wave- 
function when particles are in different states (thermal 
gas) or in the same state (BEC) (as explicited in Annexe 
I). It corresponds to the difference in the second order 
correlation function for a thermal gas and for a BEC. 

When aa > a cancelation occurs, for each channel 
of dipolar relaxation, for a given value of fc/, i.e. for a 
given value of the magnetic field. These two separate 
cancelations result in a dip in the dipolar relaxation loss 
parameter as a function of magnetic field. This dip is 
the main feature missing when calculating the cross sec- 
tions for dipolar relaxation while ignoring the effect of the 
molecular potentials, and was already discussed in [23j |. 
although never observed experimentally to our knowl- 
edge. Figure [5] shows the results of f3®*i 2 + fi^ 2 from 
eq. (IT71) and (fI5]l along with the dipolar relaxation rate 
parameter calculated from eq. © and (ITU)) . Figure [5] 
illustrates both the resonant dip in dipolar relaxation, 
and, far from the dip, the reduction by a factor of two 
in BECs due to the difference between the second order 
correlation function of a thermal gas and of a BEC at 
short range. 

To derive more physical intuition from equations (|17[) 
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the value of Fi n (r) for r 
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is proportional to 
Dipolar relaxation 



rates will therefore be proportional to the probability of 
presence of two particles at distance 
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Rdr is close to the interparticle distance R* at which the 
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FIG. 2: Comparison bet ween th e a nalytical calculation taking into 
account the effect of ae (eq. I17l and ll8t and calculations with no molec- 
ular potentials (eq. [9] and [Tot . 



molecular potential of the output channel crosses the in- 
put channel molecular potential: R* = Ti J Jj^J ^jy, (see 
Figure [I]). Despite the long range character of dipolar in- 
teractions, dipolar relaxation appears as localized, taking 
place at a specific interatomic range Rdr. oc \j\f~B. 

At very small magnetic fields, when kfag <C I, our 
results are close to the ones given in [2l[ (except for a 
factor of 2, as explained above), kfae <C 1 corresponds 
to a situation where dipolar relaxation occurs at a dis- 
tance Rdr much larger than a^. We therefore indeed ex- 
pect the molecular potentials to have no effect on dipolar 
relaxation: for I = the zero energy collision wavefunc- 
tion at distances r >> is flat and independent of a§. 
In addition, at low energy the phase shift associated to 
collisions in I = 2 vanishes. 

When the magnetic field gets larger, the effect of the 
scattering length on the incoming wave-function starts 
playing a role in dipolar relaxation. In our model, the 
cancelation of dipolar relaxation in a given channel oc- 
curs when Rdr = i.e. when dipolar relaxation occurs 
at an interatomic distance at which lies the last node in 
the incoming wave-function. However, as in the case of 
chromium a§ ~ R v dw, the analytical forms for the colli- 
sion wave functions given by eq. (|14l 1151) are not valid in 
this region. At short internuclear distances (most impor- 
tant at large magnetic fields), one needs to numerically 
calculate the incoming and output wavefunctions. We 
describe in the next part of this paper our experimental 
data and how we perform such numerical calculations to 
fit the results of such experiments. 
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FIG. 3: Atom number and size of the cloud as a function of the delay 
time t spent in ms — 3. 



MEASUREMENTS OF DIPOLAR RELAXATION 
AND SCATTERING LENGTHS 

We now turn to the experimental study of dipolar re- 
laxation. We shall see that a comparison between our 
experimental data and a numerical model taking into ac- 
count the effect of the molecular potentials provides new 
estimates of Cr scattering lengths. This method has sim- 
ilarities with those used to interpret photoassociation ex- 
periments, in which modulations in the intensity of the 
photoassociation lines reflect the structure of the collid- 
ing wavefunction [28j j , which has been used for a determi- 
nation of scattering lengths [29]. Moreover, our method 
probes the collision wavefunction at any interatomic dis- 
tance, because dipolar relaxation can be measured at any 
magnetic fields, whereas photoassociation relies on a res- 
onant condition with a molecular excited state, so that it 
only provides discrete information on the collision wave- 
function. 

To study dipolar relaxation in a ms = 3 BEC, we first 
produce Cr BECs in the ms — —3 state, as described in 
[llj . BECs are produced in a crossed optical dipole trap. 
Because trapping is weak at the end of the evaporation 
stage, we also use a magnetic field gradient along the ver- 
tical direction to compensate for gravity. After the BEC 
is produced in the crossed optical dipole trap, we recom- 
press the trap to a value such that magnetic levitation 
is not necessary anymore to compensate for gravity. We 
then remove the vertical magnetic gradient, and we ap- 
ply a bias field in the horizontal plane. We have checked 
that the number of atoms in the BEC - typically 10 000 
- is independent of the value of the final bias field. For 
this part of the paper, the value of the bias field is such 
that the Larmor frequency of the atoms is larger than 
the trap depth: dipolar relaxation results in losses. 

We then apply a pulse of radio- frequency magnetic field 
linearly polarized along the vertical direction. The rf 
frequency is swept across the rf resonance (ie around the 
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Larmor frequency). The rf power and the duration of the 
pulse (typically 2 ms) are such that more than 95 percent 
of the atoms end up in the m s = 3 state, as measured 
by Stern and Gerlach experiments. 

After a given delay time t, we apply a second rf sweep 
to the atoms, to bring them back to the ms — —3 state. 
We then abruptly (within 300 ns) release the optical trap, 
and take an absorption image after a few milliseconds 
(typically 5 ms) of expansion. If the delay time t is short 
enough, we recover a BEC, with the same number of 
atoms and the same Thomas-Fermi radius than the ini- 
tial BEC. As we do recover a BEC for delay times con- 
siderably longer than the oscillation time in the trap, our 
experiment is (to our knowledge) the first experiment to 
date to have produced a BEC of chromium atoms in the 
ms = 3 state. 

As the delay time t increases, the number of atoms 
decreases due to dipolar relaxation (see fig [3]) , without 
noticeable heating. This contrasts with what is typically 
observed in a trapped thermal gas: in that case, loss of 
atoms due to a density-dependent process mostly occurs 
at the center of the trap, so that the average energy lost 
per lost atom is less than the average energy of an atom in 
the sample. Because of this, inelastic losses in a trapped 
thermal gas results in heating of the cloud. In a BEC 
the chemical potential is flat over the size of the BEC, 
so that the energy lost per lost atom is independent of 
the position of this atom in the trap. Losses do not lead 
to heating, provided the size of the cloud adiabatically 
adapts to the change in number of particles, i.e. provided 
the loss rate is small compared to the rate of oscillation 
in the trap. The fact that we observe losses without 
heating is therefore an indirect signature of the fact that 
the chemical potential is constant over the size of the 
cloud, a pre-requisite for phase coherence in the gas. 

As the number of atoms decreases, the Thomas Fermi 
radius decreases accordingly (see the inset of fig [3]). Af- 
ter typically 20 ms, the number of atoms gets so small 
that the thermal fraction increases, and we eventually 
completely lose the BEC: this leads to an increase of the 
width of the cloud after some time of flight (see the inset 
of fig [3]) . A typical set of data is shown in fig [3] We an- 
alyze such data for delay times t short enough that the 
BEC is conserved, with no discernable thermal fraction 
(30 ms in the present case). The local decay equation 
reads: 



dn a 2 



Tin 



(20) 



where j3 r is the loss parameter associated to dipolar re- 
laxation, and Ti ps 0.1 s _1 is the one-body loss coefficient 
due to collisions with the background gas, which is inde- 
pcndantly measured. Typically, fi r n is so large that T\ 
can be neglected in eq. ([20]) . To relate this equation 
to the decay equation for the total number of atoms, we 
assume that all atoms are in the BEC, in the Thomas 
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FIG. 4: Dipolar relaxation rate parameter as a function of magnetic 
field for a ms — 3 BEC. Squares: experimental data points. Dashed 
line: calculation in the fi rst o rder Born approximation, with no molec- 
ular potential (cq. [91 and [TOII : the calculated rate parameter is divided 
by 2 to take into account the fact that the experiment is performed 
with a BEC. 



Fermi regime; we also assume that the scattering length 
in the ms — 3 state is identical to the scattering length in 
the ms = —3 state. Indeed, atoms being in the stretched 
state for both cases, the collision is purely in the 5 = 6 
molecular channel, and both states therefore have the 
same scattering length etg provided there is no Feshbach 
resonance in this range of magnetic fields. Integrating 
eq. (|2T)|) over the Thomas-Fermi density profile (and as- 
suming that there is no thermal fraction), we find that 
the equation for the total number of atoms reads: 



at 
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where a = ^f^- ( |^=J (see for example [3J|). uJ is 
the geometrical mean of the oscillation frequencies at the 
bottom of the trap, independently measured by paramet- 
ric excitation, ae = 112 ao was measured in [271 ]. There 
is an analytical solution for equation [21] which we use to 
fit our experimental data (sec fig [3] for example) in order 
to extract a value of f3 r for each magnetic field. Results 
of this procedure are presented fig [4] as a function of the 
magnetic field. As predicted in the previous part of this 
paper, we observe a dip in the experimental dipolar relax- 
ation rate parameter as a function of magnetic field, due 
to the effect of molecular potentials. As a first step to 
quantify this effect, we compare our experimental data 
to the results of the analytical model presented above. 
In the limit of zero collision energy, the analytical model 
offers direct link between and the magnetic field at 
which losses are minimum 

16(2 + S)h 



a 6 = 



3y / mgsHBB min Tr(2\/2 + S) 



(22) 
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We therefore may deduce from our measurement a first 
value of 0,6 = 117 do- However, since this value is not 
much larger than R vd w the analytical model is not valid, 
as explained above. We therefore now turn to our numer- 
ical simulations, to accurately interpret the experimental 
data. 

In our numerical simulations, the molecular potential 
is accounted for in the whole asymptotic region, and for 
both the initial and final states. The atomic spins are 
coupled, with total spin St = S\ + S2 and the dipole- 
dipole interaction can be written as ([3l|) 



H dd = 



a, 



with 



C d 



(2) 



(23) 



(24) 



where the projection of the tensor rank is taken on the 
internuclear axis. Separating angular and radial parts, 
as allowed by the Born-Oppenheimer approximation, the 
relevant initial state of an atomic pair can be written as 



\in>=\S t ,M t ,e,rni> F in {r), 



(25) 



and the j th final molecular state after dipolar relaxation 
is given by 

\out^ >= \S [ t 3 \M { t j \^\m^ > F { £ t {r), (26) 



with M t (i) + m K p = M t + mi and - l\ = 0, 2 with 
the exception of fS^ = £ = 0. For two ms = 3 atoms, 
S t = 6, M t = 6, and in the case of an initial BEC one 
additionnally has I = 0, mi = 0. 

For the corresponding final state, one has l™' = 2 
and there are two possible St values, 6 and 4. The 
initial state is thus coupled with three final molecular 
states: |oui (1) >= |6,5,2,1 >, with AE^ = gs^BB, 
and \out^ >= |6,4,2,2 >, \out^ >= |4,4,2,2 >, with 
AE (2S> = 2g$HBB. State |1 > introduced above for 
our simplified analytic model corresponds to \out^ >, 
whereas |2 > is a linear superposition of \out^ > and 
\out^ >. 

The angular part of the matrix element of eq. (|23[) is 
evaluated using standard tensor operator technique, re- 
membering that the angular momentum projections ap- 
pearing there are taken on a fixed axis (determined by 
the static magnetic field). A general formula is given in 

m. _ 

the radial wavefunctions are described in 



As in ref 31 



a simple, purely asymptotic model, based on the concept 
of nodal lines [32j, [33J . The Schrodinger equation for the 
radial collision wavefunctions is solved in the asymptotic 
part only by inward numerical integration starting from 
large r values. The constraint is that the wavefunction 
vanishes at a given nodal position, which characterizes 
entirely the inner part of the potential (in the spirit of 



the quantum defect theory) . The model depends thus on 
the asymptotic potential, i.e. essentially on the van der 
Waals Cq constant (and, to a much smaller extent, on 
the value of Cg), and on the location of the chosen nodal 
lines. It can be shown that, in a first approximation, the 
position Rq of the nodes depends both on the collision 
energy e and on the rotational angular momentum £ with 



R Q = R 0Q + Ae + Bl(£+l). 



(27) 



The constants Rqq, A and B characterize the inner part 
of the different molecular potentials, and are different for 
each St value. They can be deduced from experiment 
and we have used here values of A and B for S t = 6 and 



St = 4 derived from the data of [27j ■ A and B act here as 
small correction terms. The Cq and C% coefficients were 



also taken from 27]. The knowledge of the constant i?oo 
is equivalent to the one of the scattering length of the 
considered molecular potential and Roo is thus used as 
an adjustable parameter. 

Taking into account the confinement of the particles, 
the initial radial wavefunction Fi n (r) behaves at large 
distance like a trapped atom pair wavefunction. At short 
distance, the interaction between the two atoms cannot 
be ignored and the wavefunction is proportional to a vi- 
brational wavefunction of the St = 6 molecular poten- 
tial; its amplitude is fixed, through the normalization 
condition, by the parameters of the long-range wavefunc- 
tions. The final wavefunction F^ t (r), is, at long-range, 

an energy-normalized free wavefunction and, at short- 

(7) 

range, a molecular wavefunction of the S t potential. 
The Fermi Golden rule for a given final state (j) reads 



r (i) 



Tl 

2?r 

"a 



< in\Vdd\outW > 



= 4rl < S t ,M t ,£,mi\C d \ S$ j) , M t W , t& , m V) > 



with 



Fin(r)F^ t (r) 



r 2 dr. 



(28) 



(29) 



The dipolar relaxation rate is the sum of the Fermi 
golden rule results for dipolar relaxation towards the 
three final states. These include the relevant different 
angular weigths for the three relevant final states, which 
are in the ratio (5.400, 0.9818, 0.8182). As a consistency 
check, we observe that ratio of the coupling to \out^ > 
to the sum of the couplings to \out^ > and \out^ > is 
exactly S, equal to the ratio between dipolar relaxation 
in channels 1 and 2 defined in our analytical model (see 
eq. [131). 

The variation of the integrals with AE^ (i.e. with 
magnetic field), governs dipolar relaxation. At very small 
values of B, the radial integral mostly builds at large 
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FIG. 5: Influence of initial and final molecular potentials on dipolar 
relaxation to state |oitt* 1 *' >. The figure shows the variation of the 
square of the radial integrals I(AE^) as a function of the energy gap 
AE {1) (in G). The case where both potentials arc ignored is the dottcd- 
dashed lino. The dashed curve corresponds to the case where the initial 
potential is taken into account: the strong diminution, as compared to 
the first curve, comes from the oscillations of the vibrational wave- 
function. The final molecular potential induces an enhancement of the 
integral (solid line), which corresponds to an increase of the density 
probability close to the rotational barrier. 
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FIG. 6; Contribution of the three different channels of dipolar re- 
laxation, to states \out^' >, \out^ 2 ' > and \out^ >, and influence of 
the Cq coefficient. The figure shows the variation of the relaxation rate 
parameter of the three final states (green dashed: \out^ >; blue long- 
dashed: to \out^ >; red dot-dashed: to \out^ >), together with their 
sum (thick solid line), as a function of the magnetic field. The central 
black solid line corresponds to the central Cq value 733 ao, whereas the 
top dotted black line corresponds to 803 ao (bottom dotted black line: 
663 ao)- The value of 04 is set to 69 ao. 



interatomic distances. The value of the integral is de- 
termined by the behavior of the wavefunction at long 
distance, and the molecular potentials do not play any 
significant role. At higher magnetic fields values, on 
the other hand, the inner part of the wavefunctions con- 
tribute, and dipolar relaxation is modified by the molec- 
ular potentials, both initial and final. These general fea- 
tures, already revealed in the analysis of our analytic 
model, still hold in our numerical model. Since this mod- 
els correctly takes into account the inner part of the wave- 
function (which the analytical model failed to do) , we will 
now use it to analyze the experimental data. 

In figure [SJ we compare radial integrals describing 
dipolar relaxation to state \out^ > calculated without 
molecular potentials with similar integrals calculated tak- 
ing into account either the initial or both initial and fi- 
nal molecular potentials. (Dipolar relaxation to other 
final states is omitted here for better clarity.) In the first 
case (no molecular potentials), the initial pair wavefunc- 
tion is a radial trapped pair wavefunction (see Annexe I) 
and the final one is an energy-normalized free pair wave- 
function given by equation (Ti"5l) . In the other cases, the 
same long-range wavefunctions are smoothly connected 
to molecular wavefunctions, with St — 6 for both the ini- 
tial and final states. The C$ value (733 atomic units) is 



taken from [27 1 



It results from figure [5] that it is important to take 
into account the molecular potential both for the initial 
state and the final state. The fact that one needs to take 
into account the molecular potential in the initial state is 
also predicted by the analytical model presented above. 
However, while the analytical model predicted a zero of 
dipolar relaxation for approximately 5.5 G, the numeri- 
cal model (ignoring the molecular potential in the final 
state) leads to a zero of dipolar relaxation in this channel 
at 12 G. This disagreement confirms that the choice of 



eq. (|14[) for the incoming wavefunction is too naive, and 
that one needs to take into account oscillations of the 
wavefunction at interparticle distances smaller than ae- 

We also see from figure that the molecular potential 
cannot be ignored in the output channel either. This 
happens when the energy gain in dipolar relaxation is 
not much smaller than the height of the centrifugal I = 2 
barrier. Then, one cannot neglect the effect of the short 
range molecular potential, so that the choice of eq. [15] 
for the output channel wavefunction is also too naive. 

In figure [6[ we show the contributions of the different 
molecular channels to dipolar relaxation, as well as their 
sum. The calculation shows a dip in dipolar relaxation, 
similar to our experimental observation. As explained 
above, both the initial and the final molecular potential 
need to be taken into account for comparison between 
theory and experiment. 

In particular we find that the value of 04 has a consid- 
erable effect on dipolar relaxation for fields larger than 
1 G. We explain this rather unexpected large effect of the 
St = 4 potential on dipolar relaxation by the presence of 
a shape resonance. We have performed calculations in a 
"universal" truncated R~ 6 potential, consisting in a R~ 6 
potential limited at short range by an infinite repulsive 
wall [35j]. In the case of chromium, we expect a £ = 2 
shape resonance for an interval of 0,4 45-70 an , which in- 
cudes the 52-64 ao range recommended in [2jJ. In this 
interval, one expects a resonant increase, by tunneling 
through the rotational barrier, of the probability ampli- 
tude at short range and therefore of dipolar relaxation 
to this state. The resonance is expected close to the top 
of the barrier and the enhancement occurs at relatively 
relatively large fields. 

Hence, fortunately, as shown in Fig [71 varying 04 in- 
duces changes in dipolar relaxation only on the large 
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FIG. 7: Influence of the ag and 04 values and comparison with ex- 
periment. The figure shows the variation of the dipolar relaxation rate 
parameter as a function of the magnetic field (in Gauss). Bullets: ex- 
perimental data points, (a) Straight lines: result of the calculation for 
various values of a®: 94 ao, 98 to 107 ao with step 1 ao, 110 and 114 ao. 
The Cq coefficient is here equal to 733 ao, and — 64ao. (b) result of 
the calculation for various values of 04 from 58 ao to 70 ao with step 
1 ao. The Co coefficient is here equal to 733 ao, and a@ — 103ao. 



magnetic field side of the dip. We therefore vary the ag 
value in the range 99 ao to 107 ao, and look for the best 
agreement between theory and experiment for the low 
magnetic field side of the dip, where 04 has little effect 
on dipolar relaxation. From this first analysis, we can de- 
duce a new value for the a§ coefficient: ag = (103±4) ao- 
Setting ag to 103 ao, we can now vary 04 to find the best 
agreement between theory and experiment on the 'blue 
side' of the dip. We thus also deduce from our analysis a 
new determination of 04 = (64±4) ao- We also varied the 
Cq coefficient in the range recommended in 27J , as shown 
in Fig [SJ and we found that its effect on the reduction 
in dipolar relaxation is weaker than the effects of (Z4 and 

. The error bars that we give in our new determination 
of (J4 and ciq include variations due to the uncertainty in 
the Cq coefficient, and numerical uncertainty. 

We can now compare our new determination of o,q to 
a value strongly constrained by the position of a Fes- 
hbach resonance in d-wave collisions measured in [3lj |. 
The node position of the St — 6 wavefunction is adjusted 
to reproduce the binding energy of the first bound level 
of the ground molecular potential of Cr2, 22.85 MHz, 
which is deduced (with an error bar of 0.04 MHz) from 
the precise measurement of the Feshbach resonance in 



d-wave at 8.155 Gauss (|3l|). This measurement deter- 
mines the scattering length of the ground molecular po- 
tential S t = 6, a 6 , with a very good precision, limited 
by the precision in the knowledge of the Cg coefficient. 
Varying C§ on the whole confidence interval given in ref 
[27l ] . 663-803 atomic units, makes a§ to vary only between 
102.2 to 103.9 ao- This is to our knowledge by far the 
most precise determination of the 5 = 6 scattering length 
of Cr, and it is in very good agreement with the value 
deduced from the present analysis of dipolar relaxation 
as a function of magnetic field. 



THE INFLUENCE OF INCOMING HIGHER 
PARTIAL WAVES IN DIPOLAR RELAXATION 



In a gaussian spherical BEC, the relative wavefunction 
of any pair of particles is also gaussian and spherical, 
which results in purely s-wave scattering. Hence, despite 
the fact that dipolar interactions are long range, dipo- 
lar relaxation in a BEC can only be due to an incoming 
partial wave I = 0. For cold thermal gases, collisions for 
short-range interactions are also only in s— wave, as, ac- 
cording to the Wigner-threshold law, the collision cross- 
section in partial wave I varies as e 21 , with the collision 
energy e: a pair of particles colliding in a partial wave 
I > cannot approach each other due to the centrifugal 
barrier. The issue is different in the case of long-range 
interactions such as dipole-dipole interactions, for which 
atoms do not need to come close to each other in order to 
interact. For example non-zero partial waves contribute 
to elastic dipole-dipole scattering even at low tempera- 
tures. In this part of the paper, we discuss the case of 
dipolar relaxation, for which we will see that, at low tem- 
perature and not too low magnetic fields, contributions 
to dipolar relaxation due to collisions in non-zero partial 
waves can be neglected. We will discuss implications of 
this feature for the case of ultra-cold fermionic dipolar 
particles. 

To investigate the role of higher partial waves in dipo- 
lar relaxation, we have repeated the measurements of the 
dipolar relaxation rate parameter for thermal gases, by 
operating at different values of the optical trap depth 
at the end of the evaporation ramp, corresponding to 
different temperatures of the gas above the critical tem- 
perature. For each value of the trap depth, oscillation 
frequencies in the trap were measured by parametric ex- 
citation experiments. After the sample has been cooled 
down to its final temperature T, we set the magnetic 
field to a given value B, and send an rf sweep pulse to 
the atoms to flip the spin of all atoms to 7715 = 3 (similar 
to what was done with BECs). We then record the num- 
ber of atoms and their temperature as a function of the 
time delay after this rf sweep. To relate eq[20]to the rate 
equation for the total number of atoms, we now assume 
Boltzmann thermal equilibrium in a non-interacting gas. 
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FIG. 8: Filled squares and circles: loss parameter as a function 
of temperature for two slightly different magnetic fields close to the 
minimum of dipolar relaxation (filled circles: 3.8 G; filled squares: 3.9 
G). The lowest temperature corresponds to a pure BEC, whereas for 
other data points the gas is fully thermal. The figure shows the factor 
of 2 reduction of dipolar relaxation when BEC is reached as well as 
the fact that dipolar relaxation does not depend on temperature for 
thermal gases. 

Therefore: 

f-^-r.N (30, 

with no the peak density of the cloud. Eq. (|30|) has an 
analytical solution which we use to fit to the experimental 
data, to deduce /?. 

We have repeated this procedure for different values 
of B, and the measured thermal loss rate parameter is 
consistently a factor of roughly 2 higher than the loss rate 
parameter that we measure for similar magnetic fields for 
a pure BEC. This factor of 2 corresponds to the factor 
of two difference between the second order correlation 
function of a thermal gas and a BEC at short range. This 
was already calculated in the case of dipolar relaxation 
of Hydrogen atoms in presence of a BEC in [3(| , and we 
present our own calculation of this phenomenon in the 
first annexe of this paper. 

The results of the measurement of j3 as a function of 
temperature appear in Fig. [5] At the lowest tempera- 
ture in this figure the cloud is entirely a BEC, whereas it 
is completely thermal for the other data points. Above 
the critical temperature for BEC, we find the loss rate to 
be relatively insensitive to temperature, and about fac- 
tor of two larger than in the BEC. Outside of the dip 
in dipolar relaxation, and at sufficiently large magnetic 
fields, we do expect dipolar relaxation to be insensitive 
to temperature, as shown by eq. © and (flTT)) . We also 
show in the first theoretical annexe of this paper that this 
property remains valid in the dip, as we experimentally 
demonstrate here. 

The data shown in Fig. [8] are taken at a magnetic 
field value corresponding to the minimum in dipolar re- 
laxation shown in Fig. |4] (3.8G). At this value, the mea- 
sured value of the loss rate parameter for the thermal 
gas is a factor of 20 smaller than the loss rate parameter 



deduced from eq. ([9]) and (fTOj) . As explained above in 
the case of the BEC, both the position and the magni- 
tude of such dip are fully explained by a theory which 
only includes s — wave scattering. The similar strong 
reduction of dipolar relaxation for thermal gases at 3.8 
G therefore suggests that other incoming partial waves 
do not play any role in dipolar relaxation for these cold 
thermal gases. We have verified by numerical calcula- 
tions (shown in Fig [9]) that non zero partial waves have a 
negligible contribution to dipolar relaxation for magnetic 
fields larger than 30 mG. The calculation is described in 
Annexe I. 

One may wonder why higher partial waves do not con- 
tribute to dipolar relaxation, whereas they do contribute 
to elastic dipole-dipole interactions. For dipolar relax- 
ation at relatively high fields, the energy in the output 
channel is much higher than the initial energy. As a con- 
sequence, ^outir) oscillates much faster than ^i n (r) at 
long distances. Most of the contribution of the integral of 
^ r i7i(r) 1 i aut(T)l / ' r3 therefore comes from the region where 
^out(f) does not oscillate rapidly, i.e. from the region of 
the classical turning point of the corresponding output 
potential at the output energy, beyond the rotationnal 
barrier. For large fields, this point is at short distances. 
If the input channel is I = 0, then the input wavefunc- 
tion is constant at short distances, and non-zero. The 
integral is then also non-zero. If the input channel is a 
partial wave I > 0, ^^(r) cx r l is very small at short 
interatomic distances, and the integral is consequently 
very small. As a conclusion, only the I = partial wave 
contributes to dipolar relaxation at low collision energy 
and high magnetic fields. 

On the other hand, when the magnetic field is very 
small, ^ ut(i~) and x f , i„(r) oscillate with similar spatial 
frequency, and the integral of ^inir^outirfl/r 3 is not 
negligible, even if the initial partial wave is I > 0. Non- 
zero partial waves therefore contribute to either dipo- 
lar relaxation at low collision energy and small magnetic 
fields, or to elastic scattering due to dipole-dipole inter- 
actions. 

We have thus shown that although dipole-dipole in- 
teractions are long-range, dipolar relaxation is a purely 
s— wave phenomenon at large magnetic fields. Dipo- 
lar relaxation will therefore be suppressed for polarized 
fermions at low temperature and large magnetic fields, 
since two identical fermions collide in odd partial waves. 
Indeed, the cross section for dipolar relaxation given by 
eq. ©, (TTU)) for fermions at low temperature and high 
fields is proportional to . This means that at low tem- 
perature, the inelastic loss parameter is proportional to 
T, and also that operating at high fields is favorable 
to reduce dipolar relaxation. This feature is much en- 
couraging for two applications: first, it should be pos- 
sible to produce stable polarized Fermi degenerate gases 
with strong dipole-dipole interactions, even in a magnetic 
trap; second, mixtures of fermions in the two lowest spin 
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FIG. 9: Calculation of dipolar relaxation to state \out^ > starting 
from cither / — (red solid line) or from higher partial waves (blue dot- 
ted line: from I — 2; green dotted line: from I — 4. Dipolar relaxation 
through even larger partial waves lead to even smaller contributions) 



states should also be stable with respect to dipolar relax- 
ation: dipolar relaxation starting from a pair of atoms in 
the two lowest spin states colliding in s— wave necessarily 
results in a pair of indistinguishable fermions in an even 
partial wave, which is excluded by the Pauli principle. 

A second important consequence of the fact that dipo- 
lar relaxation is dominated by s— wave collisions at low 
temperature is that one can use the measurement of dipo- 
lar relaxation versus magnetic held B as a probe for mea- 
suring the number of atomic pairs in the gas at a typical 
distance Rdr(B). In practice, one should use numeri- 
cal simulations to take into account the inner molecular 
potentials. Measuring the dipolar relaxation loss param- 
eter as a function of magnetic fields may then be a good 
way to measure the second order correlation function in 
thermal gases as a function of interparticle distance, in 
analogy with the Hanbury-Brown and Twiss effect [36| . 
In this way, dipolar relaxation could be a probe for corre- 
lations up to distances corresponding to magnetic fields 
at which higher partial waves significantly contribute i.e 
up to about 100 nm. 



REDUCTION OF DIPOLAR RELAXATION BY 
CONFINEMENT 



It is often stated that collisions in traps are barely dif- 
ferent from collisions for homogeneous systems, and this 
is most of the times correct. However, scattering may 
be modified in certain situations of very tight confine- 
ment. Modifications of elastic scattering by confinement 
is typically obtained when the vibrational frequency of 
the trap is on the same order of magnitude or larger than 
the binding energy of the last molecular bound state in 
the collisional channel. For s— wave collisions in a har- 
monic trap, this happens when the harmonic oscillator 
size of the trap is smaller than the scattering length of 
the atoms. How confinement modifies elastic scattering 
has been theoretically investigated in the case of tight 



isotropic 3D confinement 37], ID confinement resulting 
in quasi 2D gases (38|, |39[ , and 2D confinement resulting 
in quasi- ID gases [40|. 

To date, the traps with the strongest confinement are 
created by means of optical lattices. The typical size 
of the vibrational ground state in the lattice sites is on 
the order of 50 nm, much larger than typical scattering 
lengths. A Feshbach resonance is therefore needed to 
enter the regime where scattering is modified by confine- 
ment. Confinement-induced molecules were for example 
observed for a ID Fermi gas near a Feshbach resonances 



41]. 



Inelastic scattering may also be modified in situations 
of strong confinement [38, SHI- In the limit of large 
energy release during the inelastic process (corresponding 
to large magnetic fields in the case of dipolar relaxation), 
the output channel wavefunction is not modified by con- 
finement. Inelastic losses are then modified through the 
modification of the incoming wavefunction due to con- 
finement, and a confinement-induced reduction of inelas- 
tic losses occurs when the scatterin g len gth is comparable 
to the harmonic oscillator size 3t| 42 1. 

Here, we show that inelastic scattering is also modified 
by confinement when the kinetic energy in the output 
channel is smaller than the vibrational frequency of the 
trap. In the case of dipolar relaxation, as the gain of 
kinetic energy AE is directly related to the local mag- 
netic field, it is interesting to work in strongly confining 
geometries and in weak magnetic fields. We then ob- 
serve a strong reduction of the dipolar relaxation rate 
parameter, which arises both because strong confinement 
reduces the density of states at the final energy, and be- 
cause the inelastic transitions occur at internuclear dis- 
tances R* comparable to the harmonic oscillator size of 
the trap. We emphasize that this reduction is obtained 
in a regime in which the scattering length of the atoms is 
much smaller than the harmonic oscillator size, and that 
it therefore does not rely on the proximity of a Feshbach 
resonance. 

To create very confining traps in one dimension, we 
focus a retro-reflected single frequency laser at 532 nm 
(1.5 W from a Verdi laser), with a waist of 40 /xm, on the 
Cr-BEC. As the standing wave is turned on sufficiently 
slowly (we use a 30 ms linear voltage ramp controlling 
an acousto-optical modulator), the BEC is loaded adia- 
batically in the first band of a ID optical lattice, whose 
depth is between 30 and 35 Er, where En = jtkx 7 ^ s 
the recoil energy, with A = 532 nm. The lattice beams 
propagate along the vertical axis (Oz). The vibrational 
spacing ujl/2-k in the lattice is on the order of 130 kHz, 
much larger than any other energy scale in the system 
(temperature, chemical potential); consequently, the mo- 
tion is frozen in one dimension. The size of the har- 
monic oscillator ground state wave function in the lattice 
d£ = y/h/muiL is about 35 nm. At this lattice depth, the 
tunneling time from one lattice site to the next is long 
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FIG. 10: (a) Absorption picture taken 5 ms after the optical lattice 
has been released adiabatically (band-mapping procedure), (b) (black) 
Cut of picture (a) along the axis of the lattice, as compared to (grey) a 
cross section of Bose-condcnsed atoms 5 ms after having been submitted 
to a brief pulse of the optical lattice, thus producing diffraction peaks 
at ±2hk^: the band-mapping procedure shows that all atoms lie inside 
the first Brillouin zone, showing no population in higher lattice bands, 
(c) Cut of picture (a) perpendicular to the optical lattice axis, from 
which we deduce the average temperature in the lattice. 



(400 ms) compared to the timescale of the experiment 
(up to 100 ms), so that we consider that the BEC is split 
into an array of independent 2D BECs. Indeed, time-of- 
flight images, obtained by letting the atoms expand after 
a sudden release of the optical lattice, show a gaussian 
profile along the axis of the lattice, and a Thomas-Fermi 
profile perpendicular to the lattice axis. 

After the BEC is loaded in the optical lattice, we re- 
duce the magnetic field, and set it to a low value, so that 
the Larmor frequency ujq = gsfJ-BB is chosen from be- 
low the vibrational spacing (130 kHz) to slightly below 
the lattice depth (400 kHz). For all experiments, the 
magnetic field vector is in the plane of the 2D gases, i.e. 
perpendicular to the optical lattice beams. We then flip 
the spin of the atoms with a radio- frequency sweep as de- 
scribed earlier in this paper. This produces an array of 
2D BECs in the ms = 3 state. Due to the optical lattice, 
the peak density of the BEC is raised up to 5 10 20 m~ 3 , 
noticeably higher than in the 3D BEC. 

We then observe a very rapid heating of the cloud, af- 
ter which we loose the BEC: a small fraction of the atoms 
undergoes dipolar relaxation; these hot atoms (their en- 
ergy is set by AE) remain trapped in the optical lat- 
tice, and thermalize with the other atoms, so that the 
temperature of the cloud rapidly exceeds the degeneracy 
temperature. This process happens within a timescale 
so short that we have not been able to experimentally 
analyze its dynamics. 



Before turning to the measurement of the dipolar re- 
laxation rate parameter for thermal gases in the lattice, 
we must first check that all particles remain in the lowest 
band of the lattice. In fact, for most of our experimental 
situations (when gs^BB > Twjl/2, see below), the en- 
ergy gained by one atom in a dipolar relaxation event is 
larger than the vibrational spacing in the lattice. Some 
of the atoms can therefore be excited to higher vibra- 
tionnal states in the lattice by dipolar relaxation. We 
thus operate in a regime where the time t spent in the 
ms = 3 state is small enough that most atoms do not 
undergo dipolar relaxation. We check that indeed at any 
given time the vast majority of the atoms are in the low- 
est energy band, by reducing the lattice depth in 100 
/is, which is slow compared to the timescale for band 
excitation in the lattice (10/xs), but rapid compared to 
the timescale for thermalization (on the order of 1 ms): 
this band mapping procedure [45| adiabatically transfers 
the quasi-momentum distribution in the lattice into the 
real momentum distribution. The population in the n th 
Brillouin zone corresponds to the population in the n th 
band. Typical results are represented figDJjl We do ob- 
serve that (within the signal to noise ratio) almost all 
atoms are trapped in the lowest band of the lattice. 

One ms after the rf sweep, the system therefore consists 
of an array of about forty 2D thermal clouds. The typical 
temperature is 300 nK, and the typical number of atoms 
in the central site of the lattice is 250. The 3D peak den- 
sity is tiq D = 1.5 1 20 m~ 3 , and the 2D peak density is 
tiq D = 10 13 m~ 2 , corresponding, in the central site, to a 
degeneracy parameter n^ D A 2 s» 1, where A = ^ 27T ^ nkEsT 
is the thermal de Broglie wavelength. The degeneracy 
parameter is therefore below the critical phase-space den- 
sity value for the Berezinsky-Kosterlitz-Thouless transi- 
tion 



44( | . and the 2D gas is thermal. Indeed, when the 



lattice is switched off abruptly and the atoms let free to 
expand for a few ms, the profile of the cloud reveals a 
gaussian velocity distribution perpendicular to the lat- 
tice. 

To estimate the dipolar relaxation rate coefficient, we 
first measure the heating rate in the array of 2D thermal 
clouds of ms = 3 Cr atoms. The temperature of the 
2D thermal gases is deduced by using the band mapping 
procedure explained above, and fitting the momentum 
distribution perpendicular to the lattice to extract an 
effective 2D temperature. This temperature is an average 
over the temperatures in each 2D cloud, weighted by the 
number of atoms in each 2D cloud. As shown in fig [TTJ 

(a) , the temperature of the cloud increases linearly at 
short times, whereas the number of atoms is constant. 
To deduce the heating rate, we can therefore measure 
the temperature of the cloud for two different (but short 
enough) durations spent in ms = 3, as shown in fig. [TT] 

(b) . This yields a fast measure of the heating rate in the 
lattice as a function of the magnetic field. 

To relate the measured heating rate to a rate param- 
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eter for dipolar relaxation, we also need to estimate the 
density profile of the cloud. We cannot individually ad- 
dress the different lattice sites to measure their number 
of atoms. The total number of atoms remains unchanged 
(see hgfTTI (a)), and the tunneling time (400 ms) is much 
longer than the timescale of the experiment (20 ms). We 
therefore assume that the number of atoms in each 2D 
thermal gas is identical to the number of atoms in this 
site once the BEC has been adiabatically loaded in the 
optical lattice, which we calculate, assuming the loading 
time is long enough that the chemical potential remains 
homogeneous across the sample. To estimate the density, 
we also measure the trapping frequencies by parametric 
excitation of the cloud with and without the lattice. We 
calibrate the lattice depth by applying a short pulse of 
the lattice light to the BEC and analysing the growth of 
the diffraction peaks as a function of time. 

Assuming that the density profile of the cloud along 
the axis of the optical lattice is identical to the one of the 
initial ms — — 3 BEC loaded into the optical lattice, and 
assuming a thermal profile perpendicular to the lattice, 
we can deduce the rate equation for the number of atoms 
in the ms = 3 state (N 3 ) by integrating eq. (l20l) over 
the density profile: 



dN 3 
dt 



= -P—(n d 1J ) 3 N 3 -T 1 N 3 



(31) 



where (riQ D ) 3 is the peak density of ms = 3 atoms in the 
central lattice site. This equation also assumes that most 
atoms are in the ms — 3 state. 

The small decrease of the number of ms — 3 atoms 
does not lead to losses, as AE^ is smaller than the trap 
depth. Rather, dipolar relaxation results in heating. The 
energy gain due to dipolar relaxation is only redistributed 
in two dimensions, perpendicular to the lattice axis, as 
shown by figure [TO] We assume that the trap is har- 
monic in these two directions, and we use the equiparti- 
tion theorem to relate the heating rate dT/dt to the rate 
of increase of energy due to dipolar relaxation: 



dE 
~dt 



2k B N 3 



dT 
~di 



7V2 



(n 3 D ) 3 N 3 + fhAE^/2) (32) 



These equations are valid as long as most atoms remain 
in the m s = 3 state. As AE^ = 2AE {1 \ and as both 
channels may be modified differently by confinement, we 
define a new dipolar relaxation 'heat' rate parameter 
PhZt = + 02- We then use eq. to relate the 

measured heating rate to Pheat ■ Results of this procedure 
as a function of the magnetic field are represented in fig. 

We also measure for low magnetic fields the dipolar re- 
laxation parameter for a 3D geometry, i.e. for a thermal 
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FIG. 11: Heating of the atoms in the lattice due to dipolar relax- 
ation, without losses, (a) Typical data for the number of atoms (black 
diamonds) and temperature (red bullets) of the cloud as a function of 
the time spent in ms — 3 in the lattice; the lattice depth is 35 En, the 
Larmor frequency is ujq —100 kHz. The straight line is a linear fit to 
the data used to infer the heating rate. The dotted line is a straight 
line to guide the eye (b) Temperature of the cloud for t — 0.5 ms (red 
bullets) and t — 20 ms (black diamonds), as a function of cj u - Dashed 
lines are guides for the eye. 



cloud without the lattice under otherwise similar condi- 
tions. For the weakest magnetic fields, dipolar relaxation 
results in heating of the cloud, which increases evapora- 
tion in the trap. We observe a loss rate which we there- 
fore attribute to the sole evaporation process. From the 
loss rate of atoms, and assuming that the average en- 
ergy lost per evaporated particle is (77 + l)fegT where r\ 
is the estimated ratio between the trap depth and the 
temperature [46| . and assuming equipartition of the en- 
ergy in 3D, one can therefore relate the heating rate to 



3D _ 
heat 



P1/2 + P2 using eq. (J3DJ: 



3fci 



dT 
'~dt 



{ 9 H B B)^-j=n {(i 1 /2 + p 2 ) - k B T 



dN 
~Ndi 



- v) 

(33) 
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This formula is valid as long as Ai?^^ 1 ' 2 " is smaller 
than the trap depth. For larger trap depths, we simply 
use eq. (I3T))) to determine /3 from the loss rate of particles. 
In 3D, we can directly relate the total loss rate parameter 
f3 to Pl® aV as in the limit where kf » h, P2IP1 = V2/3: 
^fa t = n(5 + 372)/?. 

We present in figure [12] (a) the experimental results 
in 3D and in 2D, 0f^ at and /3 2 ^ t , along with the result 
of the calculation for the rate parameter in 3D, corre- 
sponding to eq. © and ([TU|. and in 2D (see below). In 
this range of magnetic field strengths, our experimental 
results in 3D are in good agreement with the theory of 



21]. In addition, we observe that when the magnetic 
field is high enough, the measured rate parameter in the 
lattice is not distinguishable from the 3D case. This is in 
agreement with predictions of [391 ]: when the scattering 
length is small compared to the size of the harmonic os- 
cillator in the lattice and when the energy of the output 
channel is large compared to the vibrational energy in the 
trap, the quasi-2D loss rate parameter is not modified by 
confinement. 

When the magnetic field is lowered, the value of the ex- 
perimental 2D rate coefficient is gradually reduced com- 
pared to the 3D case. We also see in Fig. [T^lthat the 
functional form of the rate coefficient as a function of 
magnetic field is different from the 3D model. We there- 
fore identify a clear and strong effect of confinement on 
the physics of dipolar relaxation in 2D. 

To interpret our experimental results, we cannot use 
the theory previously developed in this paper. The main 
reason is that, as the magnetic field is reduced, dipolar 
relaxation takes place at longer and longer interatomic 
distances; for a magnetic fields B ss Tiojl/ gs^B, the dis- 
tance at which dipolar relaxation occurs Rdr (see eq. I19[) 
is of the order of the harmonic oscillator size in the lattice 
cll = \/ ■ As a consequence, the spherical symme- 
try of the problem is broken, and the partial wave basis 
which we have used in 3D is not a good basis anymore: 
the remaining good quantum number is that associated 
with the projection of the angular momentum along the 
axis of the ID lattice, m z . Rdr ~ a>L also corresponds 
to a situation where the energy released in a dipolar re- 
laxation event, AE^ 1 ', is of the order of Tuujl. We have 
therefore developed a 2D scattering theory based on the 
Born approximation, which includes both channels of 
dipolar relaxation, and describe explicitly the possible vi- 
brational excitation in the lattice, from the initially pop- 
ulated v — state to the v = 0, 1 and 2 states. We have 
not calculated the excitation to higher bands, because 
this would correspond to magnetic fields larger than the 
ones experimentally investigated. In this theory, lattice 
sites are assumed to be independent, and the potential in 
a given lattice site is approximated to a parabola. Our 
theoretical model is described in the second annexe of 
this paper. 
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FIG. 12: (a) Dipolar relaxation loss parameter as a function of mag- 
netic field, in the lattice (filled triangles) and without the lattice (open 
triangles). The dashed vertical lines represent (left) the magnetic field 
above which excitation by dipolar relaxation of the first vibrational 
band becomes energetically allowed (first threshold) ; (right) the sec- 
ond threshold. The red solid line is the result of the 2D theory for the 
specific orientation of the magnetic field in the experiment, and the 
black dotted line the result of the 3D theory, (b) results of the 2D 
theory for a magnetic field parallel to the 2D planes (black line) and 
perpendicular to them (grey line). 



Results of the calculation are shown in figure [T2l (b). 
The calculation shows two consecutive thresholds as the 
magnetic field is raised. Indeed, if the magnetic field 
is low enough, population of the v = 1 and v = 2 
states through dipolar relaxation is energetically forbid- 
den. When gs^sB > Twjl/2, excitation to v = 1 be- 
comes energetically allowed through channel 2 of dipolar 
relaxation: the pair of atoms gains an energy larger than 
the vibrational excitation gap of the motion in the lat- 
tice. The calculation hence shows the opening of different 
dipolar relaxation channels (v = to v = 1, then v = 
to v = 2) as the magnetic field is increased. This illus- 
trates the effect of the reduction of the density of states 
on dipolar relaxation as confinement is increased. 

As also shown in figure [12] (b), dipolar relaxation in 
2D strongly depends on the magnetic field orientation 
8 with respect to the 2D planes of the ID lattice sites 
(6 = corresponds to B parallel to the planes) . This is a 
direct consequence of the breakdown of spherical symme- 
try, and the result of an interplay between the asymmetry 
of the trap and the anisotropy of dipolar relaxation. 

In our experiment, the applied magnetic field is paral- 
lel to the 2D plane of the ID lattice sites. However, for 
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the lowest magnetic fields, our knowledge upon the ori- 
entation of the total magnetic field is reduced, because a 
stray magnetic field perpendicular to the plane may re- 
main. We therefore plot in figure [TJ] (a) , along with the 
experimental data, calculations corresponding to a mag- 
netic field consisting of a combination of a small perpen- 
dicular magnetic field and of a parallel magnetic field (the 
one we apply), which largely dominates for B > 0.025 G. 
The agreement between the experimental data and the 
theory is then good, and best for a remaining perpendic- 



ular magnetic field of B p 



0.01 G. There seems to be 



a small systematic shift between experiment and theory, 
which may be a sign that we overestimate the density in 
the lattice. 

An interesting feature in figure [T2l is that as the mag- 
netic field is reduced below Uujl/2, instead of decreasing 
as one would naively guess, the dipolar relaxation rate 
parameter increases again. As shown by our theoretical 
calculation, this phenomenon depends on the orientation 
of the magnetic field relative to the plane defined by 8. 
Indeed, in the zero magnetic field limit, dipolar relax- 
ation vanishes for 8 = tt/2, whereas it reaches a fixed 
value for 8 = 0. This is because in the first case, an 
m z = collisional wave is only coupled to m z = 2, and 
the overlap between these two waves gets smaller and 
smaller when the magnetic field is reduced, because the 
centrifugal barrier in m z = 2 strongly reduces the ampli- 
tude of the wavefunction at small interparticle distances. 
On the other hand, when 8 = 0, m z = is coupled 
to itself by dipolar relaxation, so that the coupling does 
not vanish when B goes to zero. Moreover, coupling of 
m z = to itself by dipolar relaxation increases as the 
magnetic field is lowered, because the radial parts of the 
corresponding wavefunctions get more and more similar 
as the kinetic energy in the output wave gets closer and 
closer to the initial kinetic energy. Interestingly, the value 
of 8 at which the 2D dipolar relaxation rate parameter 
is maximum at B = is different from 8 = 0. The non- 
zero value of the dipolar relaxation rate parameter as B 
goes to zero, as well as the discontinuity when the second 
threshold is reached, are both discussed in Annexe II. 

The fact that dipolar relaxation reaches a finite value 
at B = (and T = 0) is in sharp contrast to the situation 
in 3D at T = 0. This difference comes from the fact that 
even at the lowest temperatures, the lattice confinement 
creates momenta (along the lattice axis) which contribute 
to dipolar relaxation. Hence, the order of magnitude of 
the 2D dipolar relaxation rate parameter at zero mag- 
netic field is of the order of the 3D rate parameter at 
zero field for an incoming wavevector ki = 1/ot (see also 
Annexe II). This also means that, although dipolar relax- 
ation is in general smaller in 2D than in 3D, the situation 
is reversed when the magnetic field is exceedingly small 
(even smaller than what we experimentally investigate 
here). 

We have thus shown theoretically and experimentally 



a strong modification of dipolar relaxation due to con- 
finement at low magnetic fields. We observe an increase 
of dipolar relaxation when B tends to zero below the ex- 
citation gap from v = to v = 1, in agreement with 
our theoretical model. In the limit where B tends to 
0, dipolar relaxation in 2D should therefore be more effi- 
cient than in 3D for certain orientations 8 of the magnetic 
field. However, in general, dipolar relaxation is reduced 
in 2D compared to the 3D case. We observe the opening 
of new dipolar relaxation channels when the magnetic 
field is raised, corresponding to the threshold for excita- 
tion to v = 1 and v = 2. For the largest magnetic fields, 
the 2D dipolar relaxation rate parameter gets closer to 
the 3D one. We observe a crossover between a situation 
(g/isB < Tilul) for which population of the higher bands 
is completely excluded based on conservation energy ar- 
guments (dipolar relaxation happens in reduced dimen- 
sions) to a situation in which the excitation of motion 
along the third direction of space becomes possible. Then 
dipolar relaxation happens in 3D, and inelastic scatter- 
ing is not modified since the scattering length is small 
compared to the harmonic oscillator size of the trap [39| . 



The confinement-induced modification of dipolar relax- 
ation arises from the combined modification of the input 
and the output wave functions and is significant when the 
magnetic field is small enough. In contrast to all other 
predicted 
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and observed 



confinement- 



induced modifications of elastic or inelastic scattering the 
modification of inelastic scattering by confinement at low 
field does not rely on a scattering length comparable to 
the typical trap size, and can be observed without the 
use of a Feshbach resonance. This provides a general 
route to reduce dipolar losses and increase the stability 
of dipolar (spinor) gases. Our work on dipolar relax- 
ation in optical lattice will soon continue by the study 
of dipolar relaxation in 2D optical lattices, where dipolar 
relaxation should completely vanish when the magnetic 
field is smaller than the excitation energy gap to the first 
band in the 2D lattice, and parallel to the ID gases pro- 
duced by the 2D lattice. Furthermore, in a cubic 3D op- 
tical lattice, dipolar relaxation should occur only when 
the Larmor frequency exactly matches the excitation gap 
in the lattice. Working in a lattice is then a good way 
to observe the Einstein-de-Haas effect with cold atoms, 
as the excitation gap in the lattice reduces the experi- 
mental complications linked to controlling the magnetic 
field to a value very close to zero. Finally, let us mention 
the possibility, not included in our theoretical model, of 
inter-site dipolar relaxation at extremely low magnetic 
fields (i.e. when Rdr ~ A). 
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CONTROL OF DIPOLAR RELAXATION BY RF 
FIELDS 

We have shown above that dipolar relaxation can be 
controlled by means of a static magnetic field, and by 
means of confinement. We now turn to experiments 
where we show that we can also control the strength of 
the matrix element for dipolar relaxation, as well as the 
energy of the output channel, by using rf oscillating mag- 
netic fields. Controlling the energy of the output chan- 
nel is particularly appealing, as it allows for the study 
of the modification of dipolar relaxation as the process 
is changed from exoenergetic and dissipative to coher- 
ent, with no release of kinetic energy, without the need 
to control accurately a magnetic field intensity close to 
zero. As this resonant regime is reached, description in 
terms of a Fermi golden rule is expected to break down, 
and in the case of a BEC, a meanfield picture suggests the 
coherent production of structures similar to vortices, an 
effect which bears similarities with the Einstein-de-Haas 
effect [l3]. 

Our idea is to use radio-frequency (rf) magnetic fields 
to dress atoms in the ms — —3 state. Then, without 
rf, dipolar relaxation is energetically forbidden at low 
temperature. With rf, we show that dipolar relaxation 
can occur between states of different manifolds (a pro- 
cess shown in Figure H3| . The amplitude of this pro- 
cess is proportional to the dipole-dipole interaction ma- 
trix element, and depends on both the frequency and 
the Rabi frequency of the rf field. In addition, the en- 
ergy of the output channel depends on the frequency of 
the rf field. Dipolar relaxation between rf manifolds was 
to our knowledge first discussed theoretically in [47|, al- 
ready observed in our experiment as briefly discussed in 
[ijj . in a situation where the rf field was perpendicular 
to the static magnetic field. Here we study the config- 
uration where the rf field is parallel to the static field 
and we provide analytical formula for rf-assisted dipolar 
relaxation. The present work also has similarities with 
previous results in our group describing rf-association in 
the vicinity of a Fcshbach resonance as an assisted Fes- 
hbach resonance , or with earlier work with Rydberg 
atoms [5p| |. 

There is an analytical expression for the internal eigen- 
states of atoms dressed by rf, when the rf field is parallel 
to the static field, for an arbitrarily large Rabi frequency. 
As the rf field is parallel to the static field, there is no pre- 
cession, and the dressed eigenstates are products of the 
internal Zeeman substate \ms), and of a translation op- 
erator acting on the rf field. Dressed eigenstates group in 
a series of manifolds, whose spin structure is identical to 
that of undressed atoms. Inside a given manifold, the en- 
ergy of the dressed eigenstate are simply msliujQ (where 
cjo is the Larmor frequency), and consecutive manifolds 
differ by the energy of an rf photon Tiuj r f. The key point is 
the following: starting from the lowest state of energy of 
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FIG. 13: Dipolar relaxation in the presence of rf from the molecular 
physics point of view, illustrated for channel 2 (sketch) ; the difference of 
energy between initial and final channels is controlled by a combination 
of the Larmor frequency and of the rf frequency and R* is modified 
accordingly. 



a given manifold (corresponding in the case of chromium 
to ms = —3), it is now energetically possible to decay 
into a state of even lower energy, in a lower manifold. In 
fact, as the spin nature of the dressed eigenstates is not 
modified by rf, the matrix element between one magnetic 
state in a manifold and another magnetic state in another 
manifold is proportional to the matrix element between 
these two magnetic states without rf. The coefficient of 
proportionality is given by the overlap between the rf 
field components of the dressed states corresponding to 
the input and output channels. 

We consider two atoms in the \ms = —3) state. The 
dipolar relaxation reaction 



-3,-3) 



1 

71 



(|-3, -2) + |-2, -3)), 



(34) 



for which the matrix element is similar to V\ (see eq. 
(O), is energetically forbidden at low temperature. How- 
ever, if one considers this process in presence of rf pho- 
tons, it becomes energetically allowed, provided it pro- 
cedes together with the absorption of a sufficiently large 
number of rf photons AN > luq/ u) r f ■ Starting from the 
|m<j — —3) state, in presence of rf, the two different dipo- 
lar relaxation channels 1 and 2 correspond to two series 
of channels for dipolar relaxation triggered by the ab- 
sorption of N rf photons: 

The series of channels corresponding to channel 1 is: 



h3,-3, N) 



1 

7T 



-2, iV') + 1-2,-3, AT')) (35) 



AE (1) = g\i B B + (N' - N)huj rf (36) 



and for channel 2 : 



|-3, -3, AO - 
AE (2) = 2g f i B B 



-> |-2, -2, TV') 
(N' - N)hoj rf 



(37) 
(38) 



17 



We now turn to the calculation of the dipolar relax- 
ation rate parameter between these manifolds. In ab- 
sence of dipole-dipole interaction, the hamiltonian de- 
scribing two particles in presence of rf reads: 



H 



2m 



A + Hr f + gfi B B z S z 



(39) 
(40) 



where H r f = Tiuia + a + \S z {a + a + ) is the Hamiltonian 
for the radio frequency coupling with 2\y/~N = Ti£l and 
fl the Rabi frequency. The key point is that, since the 
static magnetic field is parallel to the rf field, the diag- 
onalization of Hq is exact for arbitrarily large rf power, 
and the eigenstates appear as a series of manifolds of 
dressed states [5 If : 



X,N) = T+\X,N) 



(41) 



X denotes the internal state of the pair of atoms. Tx = 
cxp (— (a — a + )) is a field translation operator, with 
Mx the spin projection of state X along the axis z. The 



eigenenergies of the dressed states 



X,N) are W x 



E x + M x huj + Nhu 



M1K. 



, where Ex is the eigen- 



energy of state X without rf or magnetic field. The term 

A/ 2 A 2 A/ 2 £"2 o 

^ — ^* is negligible compared to uj r f since for 
our typical rf power, rs U) r f, and N >> 1. 

In the spirit of the calculation presented in the first 
part of this paper, without rf, leading to eq. Eland EH we 
first neglect all molecular potentials. The dressed eigen- 
vectors are the products of three components: a plane 
wave (for the spatial wavefunction); the internal wave 
vector of each atom; the (spin dependent) translation 
operator applied to the rf field state: 

(r\k,m 1 ,m 2 ,N) = 
exp(ifc • r) exp [A(mi + m 2 )(a — a + )] \N, mi, TO2)(42) 



and the eigenvalues are 



+ NJki; - *W"'> 3 + gn B B z {mi + ™ 2 ) (43) 

Using first order perturbation theory, we now calcu- 
late the coupling matrix element between two different 
dressed states by dipole-dipole interactions: 

A = (k' ,m' l ,m' 2 ,N'\V dd \k,m 1 ^ 2 ,N) 
= (k', mi , m' 2 , N'\ TV dd T+ \k, mi ,m 2 ,N) 
= (k' ,m' 1 ,m' 2 \Vdd\k,mi,m 2 ) 

a)) \N) (44) 



x (A^'lexp —— {a 



fi'jj 



where Am = m[ + m! 2 — m\ — m 2 . The first part of the 
latter equation is identical to the one describing dipolar 
relaxation without rf: rf cannot induce coupling between 
internal states that are not otherwise coupled by dipole- 
dipole interactions. The idea of the calculation for the 
second part of this equation is to expand the exponential 
in series, then recognize the series expansion of Bcsscl 
functions when the photon number is very large. We 
thus obtain: 



(-V -M>('^(a + -»)) V) 

= Jn-n> 



f2i Am 



(45) 
(46) 



which generalizes earlier results which can be found in 
51| . The final result of our calculation is therefore simply 



related to that given in [21 



ai = ^&(?f(kf /h)kf /h\JN-N'^W\ 2 (47) 
8tt 



tr 2 



s'Wf{kf/k i )ky'/k i \j N _ N ,(2n 1 /uj)\' 2 {48) 



..(2) 
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with a = f^irV Like in the situation with no rf, fc, 
and kf are related through the conservation of energy 
equation: 

h 2 {k ( pf 



2J.2 



fifk 



+ hu>(N — N') — HBBAm (49) 



In order to obtain the final result one has to sum over 
all possible channels, i.e. over all final photon number 
states and final internal states. The result of this sum- 
mation is given, with no free parameter, in fig. 1141 along 
with experimental data. Because we measure loss param- 
eters, and not directly cross-sections, we do not directly 
plot cross section in figure [Til but the corresponding loss 
rate parameter j3 = 2aki^j 2 . As the rf frequency is 
much larger than the initial kinetic energy of the atoms, 
we also set f(kf/ki) — 2. 

To measure the rf-assisted dipolar relaxation rate co- 
efficient, we first produce either a m$ — —3 BEC with 
the same experimental procedure as above, or a thermal 
gas just above the critical temperature for BEC (by in- 
terrupting the evaporation ramp just before BEC). We 
then recompress the optical trap (as above, to compen- 
sate for gravity without applying any vertical magnetic 
field gradient), and apply a weak magnetic field, either 
vertical or horizontal. We then apply a vertically po- 
larized rf pulse, which is short compared to the lifetime 
of the BEC (a few seconds), but long enough to induce 
losses. Then we deduce the rf-assisted loss parameter by 
comparing the number of atoms with and without rf, and 
taking into account the density distribution in the trap, 
using eq. (j2"Tj) for a BEC or eq. (|3"0"]) for a thermal gas. 
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FIG. 14: rf-assistcd dipolar relaxation loss parameter as a function 
of the Rabi frequency, for a thermal gas, in Log scale. Inset (in linear 
scale): rf-assistcd dipolar relaxation loss parameter as a function of the 
Rabi frequency, for a thermal gas (black filled circles) and for a BEC 
(red squares). The measured loss rate parameter is roughly a factor of 
two smaller for a BEC than for a thermal gas. Solid line: theoretical 
results for a thermal gas. 



Typical results as a function of rf power are presented 
in fig [14j for a thermal cloud and for a BEC. For the 
data of this figure, the rf frequency is set to 300 kHz 
and the rf field is parallel to the static field. As al- 
ready observed in [48[ , the loss parameter saturates when 
Q/ui r t > 1.5. Experimental results for the thermal gas 
are compared to the predictions of eqHll ESI and we find 
a very good agreement, with no free parameter. (The rf 
Rabi frequency is independently calibrated by measur- 
ing on-resonance Rabi oscillations.) To emphasize the 
overall good agreement between theory and experiment, 
over more than two orders of magnitude in the loss pa- 
rameter, we plot both the experimental data and theo- 
retical results in Log scale for the thermal gas. As the 
inset shows, the agreement between theory and experi- 
ment is however not perfect. A disagreement of about 
20 percent is consistent with our systematic uncertain- 
ties (combined uncertainties in trapping frequencies and 
number of atoms). We also observe that the loss param- 
eter in the BEC is systematically a factor of two lower 
than in the thermal expected (see above). 

In our theoretical model, based on the Fermi golden 
rule, the loss-rate parameter is proportional to the den- 
sity of states at the output energy, as given by eq. (|36|) 
and (f38]) . As the output energy is set by the relative 
value of the Larmor frequency and the rf frequency, this 
provides a means to experimentally control it. Although 
we have not directly measured the output energy of the 




FIG. 15: rf-assistcd dipolar relaxation loss parameter for three dif- 
ferent rf frequencies (symbols: experimental data; lines: theory). 
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FIG. 16: rf-assistcd dipolar relaxation loss parameter as a function 
of Rabi frequency, for two different values of the bias field (circles: 
ujq —100 kHz ; diamonds: cjq —1.21 MHz). Losses arc much suppressed 
when the rf frequency is much smaller than ujo, indicating that rf- 
assistcd losses with the absorption of a small (< 5 in the present case) 
number of rf photons becomes energetically impossible. Solid lines : 
corresponding theoretical results 



rf-assisted losses, we have performed systematic studies 
of the rf-assisted dipolar relaxation loss parameter as a 
function of the rf power for different rf frequencies, all be- 
ing much larger than the Larmor frequency. As shown in 
figure IT51 we find a reasonnably good agreement between 
theory and experimental data, for all rf frequencies. (In 
fig [T5l the theoretical predictions are scaled down by a 
factor of two to take into account the fact that the data 
were taken with a BEC.) We hence do verify that losses 
increase as the rf frequency gets larger, due to the fact 
that the density of states in the output channel increases 
with energy. 

To further check our ability to control the energy of 
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the output channel by modifying the values of the Lar- 
mor frequency and of the rf frequency, we have performed 
experiments where the rf field is parallel to the static 
field, and the rf frequency is smaller than the Larmor 
frequency. Starting with a BEC of ms — — 3 atoms, it is 
then energetically impossible to trigger rf-assisted dipo- 
lar relaxation with the absorption of only one rf photon. 
We performed experiments with an rf frequency of 300 
kHz in a static field of .43 G, corresponding to a larmor 
frequency of uio =1.21 MHz; then, at least five rf photons 
need to be absorbed to provide enough energy for dipolar 
relaxation. We measured the rf-assisted loss parameter, 
and our results are represented in Fig [16] In this fig- 
ure, we compare the experimental measurements for a 
low (ujq =105 kHz) and a high (luq =1.21 MHz) static 
field, to the theoretical prediction summing eq. (|4T|) and 
([48]) over all energetically allowed rf-assisted dipolar re- 
laxation channels (i.e. N — N' > Amgs^BB /Tilo). 

The experiment shows that indeed losses are greatly 
suppressed in a large static field. This is because dipolar 
relaxation is energetically forbidden with the absorption 
of less than 5 rf photons, and because at low rf power, 
the Bessel functions corresponding to dipolar relaxation 
with the absorption of 5 or more photons (see eq. (|47|) 
and (|48| ) are almost equal to zero. We do observe losses 
when the rf power gets larger (Q/u > 1), and they are 
larger than what our theory predicts. We interpret this 
disagreement by the following argument: the rf field is 
not exactly parallel to the static field. Therefore, rf can 
induce spin precession; the small fraction of the atoms 
transfered in ms — — 2 can undergo dipolar relaxation 
(within their own manifold), and thus be ejected from 
the trap. 

We have thus shown that dipolar relaxation with the 
absorption of one photon can be suppressed when the rf 
photon frequency is smaller than the Larmor frequency, 
because such dipolar relaxation is then energetically for- 
bidden. This is an indirect signature that the output 
energy of a dipolar relaxation event is indeed set by the 
difference of energy between the rf frequency and the Lar- 
mor frequency. Future studies will show whether a good 
control of the output energy, by varying the rf frequency 
for a fixed bias magnetic field value, is possible or not. 
This would have the advantage of allowing the analysis of 
dipolar relaxation in a resonant regime, with the possible 
observation of the Einstein-de-Haas effect [l6| , easing the 
requirement of a very good control of the magnetic field 
at very small values. 

As a complementary information, we compare the mea- 
sured rf-assisted dipolar relaxation loss parameter when 
the rf is parallel or orthogonal to a small bias field in 
figure [TTJ In the case where the static magnetic field is 
perpendicular to the rf field, we observe an additional 
resonant loss peak (in the present case for fl/ui ~ 0.4). 
We have attributed this resonance as being due to tech- 
nical noise in the rf amplifier that we use, which produces 
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FIG. 17; rf-assistcd dipolar relaxation loss parameter for two different 
orientation of a low static magnetic field (black circles: static field 
parallel to rf field; red bullets: static field perpendicular to rf field): 
lines arc guides for the eye. 



heating and losses when a peak of the noise spectrum of 
the amplifier corresponds to the (Rabi frequency depen- 
dent, see [48]) Larmor frequency of the atoms. 

Overall, apart from the aforementioned loss peak due 
to technical noise, we see no noticeable difference whether 
the rf is parallel or perpendicular to the small bias field 
(see fig. [T7)) . As we operate in a regime where the bias 
field is much smaller than the rf field, this result is nat- 
ural: in the low bias field limit, the orientation of the 
field does not matter. However, we stress that in general 
the relative orientation of the rf field and the bias field 
does matter. In general, when the bias field is not very 
small and the rf field is perpendicular, the loss rate pa- 
rameter is different from the results of eq. (14T1) and (|48l) , 
and even the mecanism for dipolar relaxation is modi- 
fied. In particular, as for example shown in (48|, losses 
correspond to dipolar relaxation between manifolds, but 
without changing the total magnetic number of the pair 
of particles. Therefore, the energy of the output chan- 
nel is in the perpendicular case completely set by the rf 
frequency. 

Hence, we have shown that dipolar relaxation can be 
engineered by rf fields, and that the rate parameter can 
be controlled by tuning the Rabi frequency of the field. 
In the case where the rf field is parallel to the bias field, 
the energy of the output channel is set by the differ- 
ence between the rf frequency and the Larmor frequency 
(for channel 1), and by the difference between the rf fre- 
quency and twice the Larmor frequency (for channel 2). 
As it is technically very difficult to control a magnetic 
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field close to the zero value at the level required to ob- 
served the Einstein-de-Haas effect with cold atoms [lij ]. 
we suggest that one could use rf-assisted dipolar relax- 
ation for this goal: indeed, controlling both a static field 
to about 1 kHz in a given direction and a rf frequency 
to the same accuracy is a relatively easy experimental 
routine. The use of low rf power fi/w < 0.1 insures that 
dipolar relaxation with the absorption of only 1 rf photon 
dominates all other loss mechanisms, as Ji(fl/u) is then 
much higher than Ji{VL/lj) with I > 1. A very favorable 
effect is that the modification of the Larmor frequency 
due to stray transverse magnetic fields is then reduced 
by a factor Ji(Sl/w). We will work towards this goal in 
future experimental studies. 



CONCLUSION 

To summarize, we have described and experimentally 
characterized various means to control dipolar relaxation 
using external fields. Let us emphasize five major points: 

(i) We have discovered a range a magnetic fields in 
which dipolar relaxation is reduced by up to a factor of 
20 compared to previous theoretical predictions [21 1. 

(ii) We have used this effect for a new accurate deter- 
mination of the S — 6 scattering length of Cr atoms. Our 
result is in very good agreement with another new result 
presented here, obtained by the analysis of a Feshbach 
resonance in d-wave collisions [3l|. Both results reduce 



the uncertainty in ag from about 15 percent 2JJ down to 
the percent level. 

(iii) We have demonstrated that dipolar relaxation is 
purely s— wave at low temperature and high magnetic 
fields. This is in contrast to elastic dipole-dipole interac- 
tions, in which all incoming partial waves contribute. 

(iv) We have shown that dipolar relaxation can be 
modified in ID optical lattices producing independent 
2D gases. The modification by confinement of both the 
initial collision wavefunction and the final one leads to 
changes in dipolar relaxation, most significant when the 
magnetic field is small. Then dipolar relaxation depends 
on the relative orientation of the 2D gases and the mag- 
netic field, and gaps in the excitation spectrum generally 
lead to a reduction of dipolar relaxation. A reduction of 
a factor of up to 7 was observed. 

(v) Finally, we have demonstrated rf-assisted dipolar 
relaxation, which provides a way to control both the dipo- 
lar relaxation rate parameter and the energy in the out- 
put channel. 

Let us finally show how the results presented here allow 
to suggest possible future experiments. 

(i) It should be possible to use the dip in dipolar relax- 
ation for other atoms (or molecules) with large enough 
dipolar relaxation rates , for example as a new means to 
measure scattering lengths, or simply to reduce losses. 



(ii) One could use the selective cancelation of the sec- 
ond channel of dipolar relaxation and use dipolar relax- 
ation in the other channel to create a source of pairs of 
correlated atoms in state (|3,2) + |2,3)) with orbital 
momentum 1 = 2. 

(iii) Since only the s— wave collision incoming chan- 
nel contributes to dipolar relaxation at low temperature, 
dipolar relaxation should be strongly suppressed for po- 
larized ultra-cold fermions, despite the presence of elastic 
dipole-dipole interactions: dipolar Fermi gases should be 
stable in magnetic traps, as well as some trapped mix- 
tures of dipolar fermions. 

(iv) As dipolar relaxation takes place at a spe- 
cific, magnetic-field dependent, interatomic distance, one 
could use it as a probe of correlation functions in either 
thermal or quantum degenerate samples. 

(v) Strong confinement modifies and generally inhibits 
dipolar relaxation. Going to 2D or 3D optical lattices, 
dipolar relaxation should vanish below a certain mag- 
netic field threshold. In particular, by loading several 
atoms into each site of a 3D optical lattice, one should 
completely suppress dipolar relaxation, except when the 
energy released in a dipolar relaxation event exactly 
matches the energy of one excited band in the lattice. 
This proposal is promising in the prospect of observing 
with cold atoms the equivalent of the Einstein-de-Haas 
effect. 

(vi) Rf fields provide a way to control dipolar relax- 
ation, and the output energy. This fact (possibly com- 
bined with the control provided by strong confinement) 
could also be an interesting experimental avenue for the 
observation of the Einstein-de-Haas effect in a gas of 
ultra-cold atoms. 
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ANNEXE I: CALCULATION OF THE DIPOLAR 
RELAXATION RATE: CONTRIBUTION OF 
INITIAL t- WAVES WITH I AND 
TEMPERATURE-DEPENDENCE 

This theoretical annexe describes dipolar relaxation 
from the many-body point of view. We use this approach 
to interpret three experimental results detailed in the pa- 
per: the dipolar relaxation rate parameter is twice higher 
in a thermal gas than in a BEC; dipolar relaxation does 
not depend on temperature at large magnetic field; at 
large magnetic fields, only s— wave contributes to dipo- 
lar relaxation. 
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We describe the initial state of the N indistinguishable 
bosonic atoms (ms = 3 Cr atoms) by a fully symmetric 
state in the occupation number representation: 



\In >= 



>= 



n 



\vac) 



(50) 



where the ttj's represent the different vibrational states 
in the trap. a' a . is the particle creation operator relative 
to the single particle state a^. \vac) is the vacuum state. 
We consider two simple situations: a pure BEC, with 
n ao = N and all other occupation numbers equal to zero, 
and a pure thermal cloud in an isotropic harmonic trap 
of frequency u>/ (2ir); the occupation numbers then follow 
the Maxwell-Boltzmann statistics, with 



N exp(-/3E ai )/ exp(-/3£ Q , ) 



(51) 



where /? = 1//cbT, ks being the Boltzmann constant and 
T the temperature. In the second-quantization formal- 
ism, the dipole-dipole interaction operator reads 

1 



V 



- 4 a l axa < < &\Vdd\K >, (52) 



where the sum runs over all atomic trapped and un- 
trapped states and where the matrix element is a non- 
symmetrized two-atom matrix element. 

In the final state of a dipolar relaxation event, two non- 
zero occupation numbers of the states oti (i — n, m) which 
are initially populated (trapped states) will decrease by 
one (or one occupation number will decrease by two), 
whereas two new (untrapped) states 7^ (i = (p,q)) will 



now have 



1 (or one new state will have n = 2) 



The final occupation numbers are thus n' a . and n^. , with 
(we give here only the case n ^ m and p ^ q, the other 
cases being readily obtained) 

n ' a — n ai f° r i ^ n and i =/= m 
n a = n ai — 1 for i = n or i = m 

n! li = 1 for i = p and i = q. (53) 

The number of initially non-populated untrapped states 
7i is very large, so that in both cases (BEC or thermal 
cloud), the terms with p — q contibute in a negligible 
way. In addition, for the thermal case, as the number of 
populated trapped states is large, the terms with n = m 
can also be safely neglected. For the thermal cloud, the 
dipole-dipole matrix element corresponding to transition 
from states a n ,a m to states 7 P ,7 9 thus reads: 



N 



exp(-P(E ari +E am )/2) 



< Out\V\In >= 
< 7 P lq\Vdd\a n ®m > (54) 



Eiexp(-^ Q J 
and, for a condensate in the state ao, 
N 

< Out\V\In >= -j= < -f p jq\V dd \a ao > ■ (55) 



In both cases, the bras and kets of the matrix elements 
are now symmetrized two-atom wavef unctions. In order 
to conveniently deal with the dipole-dipole interaction, 
we introduce the coupling of the spins, with total spin 
St = S1+S2, and the separation of the center of mass and 
relative motions. Instead of the two-atom wavefunctions 
of equations (|54|) and (|55|) . which are functions of the 
individual atom positions R\ and R\, we use atom pair 
wavefunctions of the new variables Rc — {R\+R<i)/2 and 
r = R\—R2- These variables are indeed separable for the 
trapped motion in the harmonic trap as well as for the 
free motion. The two-body states, with wavefunctions 
4>c and (f> re i describing the center-of-mass and relative 
motions, can be formally written in terms of the old ones 

|TOsi,ai,m S 2,Q!2 >= \St,M U (j}cAr<il> 

<Pc,<Pr<,l 

< S t ,M t ,cj)C:4> r ei\m S i,ai : ms2,a2 >, 

with E^ c + E^ rel = E ai + E a2 (56) 

The dipole-dipole matrix element is given by equations 
very similar to ([54]) and ((55]): 



< Out\V\In >= N 



exp(-/3(^+^ | )/2) 



< S' t ,M{, ^^flVddlSuMt, <j?S4>7ei >, (57) 
for the thermal cloud, and 

N 



< Out\V\In >-- 



V2 



x < Si, Mi <j>°c€tl\Vdd\St, M t ,C€ >, (58) 

for the condensate. The ket in the latter equation is 
readily obtained, by variable change, from the two-atom 
wavefunction lao^o >• I n both cases, the dipole-dipole 
interaction does not concern the center of gravity mo- 
tion and one has 0™* = 0™. The decay rate, through 
dipolar relaxation, of the initial N-atom state is given by 
the Fermi Golden Rule (following first order perturbation 
theory) 



2tt 



r= — J2\< Out\V\In>\ 



(59) 



(60) 



and is related to the atom loss by 

r = --— 

2 dt 

In the BEC case one obtains 

X^-E l< S' u M' t ^° r lt\V dd \S u M u ^ el >| 2 (61) 



r 



where out in the sum stands for S' t , M' t , 4>Vei- Assuming 
that dipolar relaxation takes place in an otherwise non- 
interacting BEC, the relaxation rate parameter is finally 
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given by 



HO ^ 

O IL I, 



Pr 



\<S' t , M[, <t>°tt\V M \S t , M t , 4>Z > I (.62) 



with clho — y^JJmuj) the characteristic length of the 
trap. The initial and final wavefunctions 4> re i contain a 
spherical harmonic, with I — mi = 0, for the initial state 
and with I = 2 and mi = M t — M[, for the final one, mul- 
tiplied by a radial function. As explained in the main 
part of the text, the radial part of the initial (resp. final) 
wavefunction behaves, at short distance, as a molecular 
vibrational function of the St (resp. S' t ) molecular po- 
tential curve. At long distance, it behaves either as a 
gaussian wavefunction for the initial state, or as a free 
pair wavefunction, for the final one (eq. 1151) . 

The case of an interacting BEC is harder to solve, as 
the relative and the center-of-mass motions of a pair of 
particles are not separable: the wavefunction describing 
the relative motion of the atom pair is not independent 
of the relative directions of Rc and r. For a BEC in the 
Thomas Fermi regime with radius Rtf, the typical colli- 
sion energy Ti 2 jmR\ F is very small; we therefore assume 
an initial I = mi = partial wave for any collision. We 
use an 'average' radial pair wavefunction which we calcu- 
late by evaluating the radial density of probability i.e. by 
the integration of the square of the two-atom Thomas- 
Fermi wavefunction over all the variables other than \r\. 
The final Thomas-Fermi pair wavefunction is taken as 
the square root of this density of probability: 



el.TF 



(R) 
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16V7 



(2 - R) V32 + 64i? + 24R 2 + 3R 3 . (63) 



with R = t/Rtf- The relaxation rate parameter in the 
Thomas Fermi regime is given by 

28^ 2 



Pr = 



-R 



3 \ ^ c/ ti j f A ou t 



out 



< S' t ,Ml,r r u J\V dd \S t ,M t ,rt n e i > 2 (M) 



where is connected at short distance to a molecular 
wavefunction of S t , and at long distance to the Thomas 
Fermi pair wavefunction (eq. I?53")) . We stress that the 
choice of the long distance radial wavefunction for the 
calculation of /3 is of little importance, as P involves an 
exact normalization to the atomic density. Hence, eq. 
([62"]) and eq. (|64|) lead to the same numerical result. 

In the case of a thermal cloud, both eq. (|59|) and (|60|) 
apply. However, the dipole-dipole interaction operator 
written in the second quantization form (eq. (|52l) ). when 
applied to the initial many-body state given by |50[ intro- 
duces another summation over the inital two-body states. 
One has 



exp(-PE in ) 



T=^N 2 V 



\<S' t ,M' t , t 



Cf\V dd \S u M u 



"rel 



> 



(65) 



where £^„ is the total energy of the initial state (in- 
cluding the motion of the center of gravity). The cen- 
ter of mass motion and the relative one are both har- 
monic, with frequency lo/(2tt). The initial pair states 
are characterized by quantum numbers Ni, Li, Mi, with 
ENi.Li — (2Ni + Li + 3/2) Tluj, for the global motion and 



,£i,mi, with E„ 



(2rii + li +3/2) Tiuj for the rela- 



tive one [52|. It is important to notice that, because of 
symmetry properties, one has to restrict to even values of 
£i, but that both even and odd values of Li are allowed. 
The sum in the denominator of equation [BS] still runs on 
atomic harmonic states. Provided that £ = {Tiuj) / '(ksT) 
is small, one has 



(66) 



The equation (|65p can be written as 
2irN 2 



r = 



hD 2 



E 



S' t , M' t ,1' ,Ni Xi , Mi ,m ,m; 

< S' t ,Ml,E',t!,m'\V M \St,M t ,n i ,li,m i > 2 



(67) 



where the energy E' of the final state is fixed by the 
resonance condition and m' = M t — M' t . The sum over 
Ni,Li, Mi of exp( — C-^Ni^Li ) is equal to D too. The atom 
loss for a gaussian atomic cloud is related to the relax- 
ation rate by 



dN 
~dt 



1 



p r N a-* C 



(4tt) 3 / 2 

The relaxation rate is thus given by 

a ( 4 ^) 5/2 



(68) 



„3 A3/2 
a HO S 



E 



exp(-C(2ra; + ^)) 



S'M'l',rii 



|< S , t ,Mi,E , ,e,m'\V id \St,M t ,n i ,li,mi >f . (69) 

The matrix element in this equation can be split into 
angular and radial parts. The angular part yields the 
selection rules: triangular relations for the sets (St, S' t , 2), 
(li, I' , 2), and M[ + m! = M t + m^. The radial integral is 



I(E') = / 
Jo 



r 2 dr. 



(70) 



As explained in the main part of the text, the radial part 
of the initial (resp. final) wavefunction are connected at 
short distance, to a molecular vibrational function of the 
St (resp. S't) molecular potential curve. 

We can now estimate the contribution of the different 
partial waves to dipolar relaxation. Figure |H] shows the 
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variation with AE of the square of the radial integrals 
for various initial and final quantum numbers. All in- 
tegrals with £i and < 10 vanish for a magnetic 
field larger than approximately 30 mG. For large enough 
energy gaps between initial and final states, the initial 
pair states with li ^ do not contribute to relaxation 
because the effect of the rotational barrier on the final 
wave-function strongly limits the efficiency of the dipolc- 
dipole interaction. As already mentioned in the main 
part of the text, although dipole-dipole interactions are 
long-range, only s-wave collisions contribute to dipolar 
relaxation at large enough magnetic fields. 

We can also use eq. ([69]) to study the temperature de- 
pendence of dipolar relaxation. It is useful to introduce 
a radial wave function of a slightly different form (corre- 
sponding to a different normalization), G(r) — rF(r), 
where F(r) corresponds to the previously used radial 
wavefunction. Then, for I = initial states, the ma- 
trix elements for dipolar relaxation are with a very good 
precision proportional to the derivative at r = of the 
long-range wavefunction corresponding to G{r), a, i.e. 
its slope. This property arises from the fact that the am- 
plitude of the wavefunction in the region where dipolar 
relaxation occurs is set by the long-distance wavefunc- 
tion describing the motion in the trap (both are con- 
nected to describe both the molecular properties and 
the motion in the trap), which is proportional to a. 
This property stands both for harmonic oscillator and 
Thomas-Fermi wavefunctions. The slope at the ori- 
gin of the harmonic oscillator can easily be evaluated: 

<JHo{ni) ~ a/| K + l/2) 1 /4 (a HO )- 3/2 (for large m), 



and a H o(0) = U 



2\ 1/4 -3/2 



l HO 



In case of a Thomas-Fermi 



pair wavefunction, <ttf = y^f R-tf^ 2 

Introducing in equation (|69|) the proportionality to 
{iii + 1/2) 1 / 4 and to (ano) 3 ^ 2 of the matrix element, 
and summing over rii, ti and rrii, we find that the relax- 
ation rate of a thermal cloud is finally given by 

Si, Ml 



< S[,M' t , E>, 2, Mt - M' t \V dd \S u M t ,0, 0, > ; 
cjhoW 



in agreement with our experimental observations shown 
in figure [8j which shows that the dipolar relaxation rate 
parameter is to within experimental uncertainty identical 
at 0.8/iK and 4.5^K. 

In figure[Sl we also experimentally verify that the dipo- 
lar relaxation rate parameter is twice higher in a thermal 
gas than in a BEC. Our theoretical treatment also ex- 
plains this fact. Indeed, the relaxation rate parameter of 
the non-interacting BEC (eq. 162]) can be written as: 



T £ 4lT 



Si Mi 



< S' t , Mj, E', 2, M t ~ M' t \V dd \S u M t ,0, 0, > 2 , 

Comparing eq. (|TTj) and eq. ([?2"]). we deduce that the 
relaxation rate parameter for a thermal cloud is twice 
larger than for a condensate, as already known from, for 
example, [3(3 ] . Such result holds in the dipolar relaxation 
dip (as shown in fig. [5]). 

For low magnetic fields, the result would be different 
for two reasons. First, partial waves I ^ contribute to 
dipolar relaxation in a thermal gas, and not in a BEC. 
Second, if the distance at which dipolar relaxation oc- 
curs is larger than the de Broglie wavelength, decoherence 
should be taken into account, and the coherent state cho- 
sen for the many-body input state (eq. ([50)) ) would not 
be a proper choice anymore. 

The relaxation rate parameter of an interacting BEC 
in the Thomas Fermi regime can also be written in a 
similar way: 



£ 4 » 

Si Mi 



< SI, Mj, E', 2, M t - M' t \V 12 \St, M t ,<hF > S 

2 



(73) 



(71) 



As the molecular wavefunction in the non interacting 
case and in the Thomas Fermi regime are proportional 
to respectively <7f/o(0) and utf (the remaining of the 
wavefunction being identical in both cases, only set by 
| St, Aft >), we also confirm, by comparing eq. (f72[) and 
eq. ([73)) that the result of the calculation of the dipolar 
relaxation rate parameter does not depend on the exact 
shape of the radial wavefunction, as stated above. 



the fraction depending only on S' t ,M' t ,E' and on the 
molecular parameters. It is now clear that the rate is in- 
dependent of the temperature of the thermal cloud, pro- 
vided E' >> E. We therefore understand why dipolar 
relaxation is temperature independent for large enough 
magnetic fields. This is in agreement with the results 
found in the framework of the first order Born approxima- 
tion, neglecting all molecular potentials [21] . This result 
remains valid even when molecular potentials are taken 
into account, and in particular in the vicinity of the dip in 
dipolar relaxation described in this paper. Such result is 



ANNEXE II: DIPOLAR RELAXATION IN 2D 

The Hamiltonian in the center of mass coordinate of 
two atoms reads 

U -HL A + Vcxt{z) + Vdd ^ + 9^B-S (74) 
m n 



where V cxt (z) = (m/4)w| / z 2 , B is the magnetic field 
S = Si + S2 is the total spin and V dd (f) is the dipole- 
dipole interaction (see eq. ([!}). wl/27t corresponds to 
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the vibrational frequency in a lattice site (see main part 
of the paper). 

The eigenstates of the Hamiltonian in absence of inter- 
action read 



(75) 



with the energy 



fi 2 k 2 
_ n k 



HujLn — g[iBB[i + 8 (76) 



where p = (x,y), k = (k x ,k y ) is the momentum in the 
xy— plane, 4> n (z) are the eigenstates of the harmonic 
oscillator in the z direction and \fx) are the eigenstates 
of the spin part of the Hamiltonian (defined in the main 
body of this paper), with \i = 0, 1, 2, and 8 is a constant 
insuring that —gfigBfi + 8 is the total magnetic energy 
of the atom pair. 

The scattering cross section in two dimensions can 
be calculated in the first order Born approximation and 
reads 



Tinlnf (&»>*/) — ( 



2nk, 



\Ei = E, 



(77) 

ki,rii, fXi and kf,nf, fif are the initial and final quantum 
numbers, and the energy must be conserved. The result 
is analytical. The previous result is correct for distin- 
guishable particles; in case of identical particles we have 
to replace 



y>¥ \VddW )±(i>v \v dd \r\ ) 

kf,fj,f ki.j-Li kf,fj,f —hi, fii 



(78) 



where the sign + refers to bosons and — to fermions. In 
order to calculate the total cross cross-section we inte- 
grate and sum over all possible final states and we per- 
form an average over the initial states with the same 
energy: 



(ki,k f ) = — 



1 



2tt 



cr^ 1 % s s (ki,kf)d9fd9i\ Ez =E f 
(79) 

(cos(9ij),sm(9ij)). 

The integral in eq. (|7T[) can be rewritten in the follow- 
ing form 



where kij = 



: > J 



(ipf\V dd \ip t ) = ^-(^f\Vd d (ki - kf,q z )\Hi) f ni ,n f (-qz)dq z 

(80) 

where V dd (q) is the Fourier transform of the dipole-dipole 
potential V dd (f) and f ni ,n f (Qz) is the Fourier transform 

Of fm,nf(z) = {<f>nf{z))*<t>ni(z)- 

An important issue is the orientation of the magnetic 
field with respect to the 2D plane. Let us consider a 
magnetic field of the form 

B = B(cos 9u x + sin 9u z ) (81) 



In this case the matrix elements read 



(MV dd (q)\0) 



(82) 



= 47r^2 (Qx sm — 1z cos 9 + iq y ) (q x cos 9 + q z sin ( 



(2\V dd (q)\0) 
sin 9 — q z cos i 



(83) 



+ iq v f 



In the scattering process the total energy is conserved and 
all processes are not allowed; in particular the following 
condition must be fulfilled 



tfk) 



2 U2 



h z k 



+ hoj L {ni - rif) - gsHBB(ni - (if) > 0, 

lit lit 

which implies that there is magnetic field dependent up- 
per bound for the final quantum number nt. In our ex- 
periment the following conditions are fulfilled: h 2 kf /m < 
Tilol and = 0. As a consequence if the magnetic field 
is not sufficiently intense the transition to higher state in 
the z-direction is forbidden. Since the difference of spin 
state can be ^ — pf = —1,-2 the first transition occurs 
for rbf = 1 and [ii — fif — —2 at a critical magnetic field 



2 U2 



2g s HBB crl = huj L - 



h k. 



(85) 



The second threshold appears for rif — 2 and fj,i — fit 
—2 at a critical magnetic fields 



2 u2 



(86) 



The next threshold appears for rif — 1 and /Uj — jif = —1 
at a critical magnetic field 



2 1,2 



gsHBB cr3 = hu) L 



h z k. 



m 



(87) 



The second and the third thresholds coincide in the limit 

— * > 0. The first threshold corresponds to a pair of 

particle gaining a kinetic energy Tuol through dipolar re- 
laxation in channel 2; after dipolar relaxation, each atom 
of this pair is then in a coherent superposition of the first 
two vibrationnal level of a lattice site. 

An important aspect that strongly affects the inelastic 
scattering is the orientation of the magnetic field 9. We 
can distinguish two limiting cases: the magnetic field or- 
thogonal to the xy~ plane, 9 = ir/2 and the magnetic field 
lying in the xy— plane, 9 — 0. In fact different orienta- 
tions of the magnetic field introduce different functional 
behaviors of cr(B), and different selection rules since cer- 
tain transitions are inhibited. The simplest example is 
the transition to rif — 1 with fii — \if = —2 for 9 = tt/2; 
in this case the matrix element (tpf\Vdd\ipi) vanishes due 
to symmetry since the integration along the z direction 
of an odd function vanishes. This is the reason why the 
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first threshold disappears in hg [12] (b) (grey curve) . An- 
other important aspect is the limit of B — > B cr : when 
9 = 7r/2, the cross section of all channels opening above 
B — > B cr vanish, and this is the reason why the cross- 
section (j(B) is then continuous. If 9 ^ ir/2 the cross 
section of channels opening above B — > B cr do not gen- 
erally vanish when B — » B cr ; as a consequence cr(B) can 
present a discontinuous behavior. In practice the aspect 
of the first and second threshold is continuous, and the 
third is discontinuous. For the limit B — > the behavior 
also differs as a function of 9. For 9 — tt/2, the dipolar 
relaxation rate parameter vanishes at low magnetic field 
and temperature T = 0, whereas for 9 = 0, it is, to within 
a numerical factor, 



[8] 
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P2D,T=0 



S 2 d 4 m 1 
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